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SUMMARY 


An interactive FORTRAN program that provides the user with a powerful and flex- 
ible tool for the linearization of aircraft aerodynamic models is documented in this 
report. The program LINEAR numerically determines a linear system model using non- 
linear equations of motion and a user-supplied linear or nonlinear aerodynamic model. 
The nonlinear equations of motion used are six-degree-of-freedom equations with sta- 
tionary atmosphere and flat, nonrotating earth assumptions. The system model deter- 
mined by LINEAR consists of matrices for both the state and observation equations. 
The program has been designed to allow easy selection and definition of the state, 
control, and observation variables to be used in a particular model. 


INTRODUCTION 


The program LINEAR described in this report was developed at the Dryden Flight 
Research Facility of the NASA Ames Research Center to provide a standard, docu- 
mented, and verified tool to derive linear models for aircraft stability analysis 
and control law design. This development was undertaken to address the need for 
the aircraft specific linearization programs common in the aerospace industry. 

Also, the lack of available documented linearization programs provided a strong 
motivation for the development of LINEAR; in fact, the only available documented 
linearization program that was found in an extensive literature search of the field 
is that of Kalviste (1980). 

Linear system models of aircraft dynamics and sensors are an essential part of 
both vehicle stability analysis and control law design. These models define the 
aircraft system in the neighborhood of an analysis point and are determined by the 
linearization of the nonlinear equations defining vehicle dynamics and sensors. 

This report describes a FORTRAN program that provides the user with a powerful 
and flexible tool for the linearization of aircraft models. LINEAR is a program 
with well-defined and generalized interfaces to aerodynamic and engine models and 
is designed to address a wide range of problems without the requirement of program 
tnodif ication . 

The program LINEAR numerically determines a linear system model using nonlinear 
equations of motion and a user-supplied linear or nonlinear aerodynamic model. 
LINEAR is. also capable of extracting both linearized engine effects, such as net 
thrust, torque, and gyroscopic effects and including these effects in the linear 
system model. The point at which this linear system model is defined is determined 
sither by completely specifying the state and control variables, or by specifying 
an analysis point on a trajectory and directing the program to determine the con- 
trol variables and remaining state variables. 

The system model determined by LINEAR consists of matrices for both the state 
and observation equations. The program has been designed to provide easy selection 
Df the state, control, and observation variables to be used in a particular model. 
Thus, the order of the system model is completely under user control. Further, the 
program provides the flexibility of allowing alternate formulations of both the 
state and observation equations. 
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LINEAR has several features that make it unique among the linearization programs 
common in the aerospace industry. The most significant of these features is flexi- 
bility. By generalizing the surface definitions and making no assumptions of sym- 
metric mass distributions, the program can be applied to any aircraft in any phase 
of flight except hover. The unique trimming capability, provided by means of a user 
supplied subroutine, allows unlimited possibilities of trimming strategies and sur- 
face scheduling, which are particularly important for asymmetrical vehicles and 
aircraft having multiple surfaces affecting a single axis. The formulation of the 
equations of motion permits the inclusion of thrust-vectoring effects. The ability 
to select variables for the state, control, and observation vectors of the linear 
models, without program modification, combined with the vast array of observation 
variables available, allows any analysis problem to be attacked with ease. 

The interactive LINEAR program is an extension of the batch LINEAR program 
described in Duke and others (1987). Data describing the aircraft and test case 
are input to the program through a terminal or formatted data files similar to the 
input files for the batch program. All data can be modified interactively from case 
to case. The aerodynamic model can be defined in two ways: a set of nondimensional 

stability and control derivatives for the flight point of interest, or a full non- 
linear aerodynamic model as used in simulations at Ames-Dryden. 


NOMENCLATURE 


The units associated with the following quantities are expressed in the English 
system. LINEAR will work equally well with any consistent set of units with two 
notable exceptions: the printed output and the atmospheric model. Both the printec 

output and the atmospheric model assume English units. 

Where applicable, quantities are defined with respect to the body axis system. 


Variables 


A state matrix of the state equation x = Ax + Bu; or, axial force (lb) 

• 

A' state matrix of the state equation Cx = A'x + B'u 

a speed of sound in air (ft/sec) 

normal acceleration not at center of gravity (g) 
a nx x body axis accelerometer at center of gravity (g) 

a nx, i x t> od y axis accelerometer not at center of gravity (g) 

a n y y body axis accelerometer at center of gravity (g) 

a ny, i y body axis accelerometer not at center of gravity (g) 

a nz z body axis accelerometer at center of gravity (g) 
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' a nz,i z body axis accelerometer not at center of gravity (g) 

a x acceleration along the x body axis (g) 

ay acceleration along the y body axis (g) 

a z acceleration along the z body axis (g) 

• 

B control matrix of the state equation x = Ax + Bu 

j • 

B 1 control matrix of the state equation Cx = A'x + B'u 

b wingspan (ft) 

C C matrix of the state equation Cx = A'x + B'u; or, force or moment 

coefficient 

c mean aerodynamic chord (ft) 

D dynamic interaction matrix for the state equation x = Ax + Bu + Dv; or, 

drag force (lb) 

D' dynamic interaction matrix for the state equation Cx = A'x + B'u + D'v 

det determinant 

E dynamic interaction matrix for the observation equation y = Hx + Fu + Ev 

E' dynamic interaction matrix for the observation equation 

y = H'x + Gx + F'u + E'v 

E s specific energy (ft) 

F feedforward matrix of the observation equation y = Hx + Fu 

F* feedforward matrix of the observation equation y=H'x+Gx+F'u 

I 

fpa flight path acceleration (g) 

G G matrix of the observation equation y = H'x + Gx + F'u 

g acceleration due to gravity (ft/sec^) 

h altitude (ft) 

H observation matrix of the observation equation y = Hx + Fu 
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observation matrix of the observation equation y=H'x+Gx+F'u 

angular momentum of engine rotor ( slug-f t 2 /sec ) 

aircraft inertia tensor (slug-ft 2 ) 

rotational inertia of the engine (slug-ft 2 ) 

x body axis moment of inertia (slug-ft 2 ) 

x-y body axis product of inertia (slug-ft 2 ) 

x-z body axis product of inertia (slug-ft 2 ) 

y body axis moment of inertia ( slug-f t^) 

y-z body axis product of inertia (slug-ft 2 ) 

z body axis moment of inertia (slug-ft 2 ) 

total body axis rolling moment (ft-lb) 

generalized length (ft) 

Mach number; or, total body axis pitching moment (ft-lb) 
aircraft total mass (slug) 

normal force (lb); or, total body axis yawing moment (ft-lb) 
load factor 

specific power (ft/sec) 
roll rate (rad/sec) 
ambient pressure (lb/ft 2 ) 
total pressure (lb/ft^) 
pitch rate (rad/sec) 
dynamic pressure (lb/ft^) 
impact pressure (lb/ft 2 ) 


yaw rate (rad/sec) 



Reynolds number 


Reynolds number per unit length (ft -1 ) 
wing planform area (ft 2 ) 

ambient temperature (°R); or, total angular momentum ( slug-ft 2 /sec 2 ) 

total temperature (°R) 

velocity in x-axis direction (ft/sec) 

control vector 

total velocity (ft/sec) 

velocity in y-axis direction (ft/sec) 

dynamic interaction vector 

calibrated velocity (ft/sec) 

equivalent velocity (ft/sec) 
vehicle weight (lb) 

velocity in z-axis direction (ft/sec) 
total force along the x body axis (lb) 
state vector 

thrust along the x body axis (lb) 

sideforce (lb) 

observation vector 

thrust along the y body axis (lb) 

total force along the z body axis (lb) 

thrust along the z body axis (lb) 

angle of attack (rad) 
angle of sideslip (rad) 
flight path angle (rad) 

displacement of engine from center of gravity along x body axis (ft) 



displacement of engine from center of gravity along y body axis (ft) 

displacement of engine from center of gravity along z body axis (ft) 

lateral trim parameter 

longitudinal trim parameter 

incremental rolling moment (ft-lb) 

incremental pitching moment (ft-lb) 

incremental yawing moment (ft-lb) 

directional trim parameter 

thrust trim parameter 

incremental x body axis force (lb) 

incremental y body axis force (lb) 

incremental z body axis force (lb) 

pitch angle (rad) 

density of air (slug/ft^) 
torque from engines (ft-lb) 
gyroscopic torque from engines (ft-lb) 
roll angle (rad) 

tilt angle of acceleration normal to the flightpath from the vertical 
plane (rad) 

heading angle (rad) 

engine angular velocity (rad/sec) 

Superscripts 

derivative with respect to time 
transpose of a vector or matrix 

Subscripts 

aerodynamic reference point 


total drag 



altitude 


L total lift 

1 rolling moment 

M Mach number 

m pitching moment 

max maximum 

min minimum 

n yawing moment 

p roll rate 

q pitch rate 

r yaw rate 

s stability axis 

x along the x body axis 

y along the y body axis 

Y sideforce 

z along the z body axis 

o standard day, sea level conditions; or, along the reference trajectory 

FORTRAN Variables 

AIX inertia about the x body axis 

AIXE engine inertia 

AIXY inertia coupling between the x and y body axes 

AIXZ inertia coupling between the x and z body axes 

AIY inertia about the y body axis 

AIYZ inertia coupling between the y and z body axes 

AIZ inertia about the z body axis 

ALP angle of attack 

ALPDOT time rate of change of angle of attack 
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AMCH 

AMSENG 

AMSS 

B 

BTA 

BTADOT 

CBAR 

CD 

CDA 

CDDE 

CDO 

CDSB 

CL 

CLB 

CLDA 

CLDR 

CLDT 

CLFT 

CLFTA 

CLFT AD 

CLFTDE 

CLFTO 

CLFTQ 

CLFTSB 

CLP 

CLR 

CM 


Mach number 

total rotor mass of the engine 

aircraft mass 

wingspan 

angle of sideslip 

time rate of change of angle of sideslip 
mean aerodynamic chord 
total coefficient of drag 

coefficient of drag due to angle of attack 

coefficient of drag due to symmetric elevator deflection 

drag coefficient at zero angle of attack 

coefficient of drag due to speed brake deflection 

total coefficient of rolling moment 

coefficient of rolling moment due to beta 

coefficient of rolling moment due to aileron deflection 

coefficient of rolling moment due to rudder deflection 

coefficient of rolling moment due to differential elevator deflection 

total coefficient of lift 

coefficient of lift due to angle of attack 
coefficient of lift due to angle-of-attack rate 
coefficient of lift due to symmetric elevator deflection 
lift coefficient at zero angle of attack 
coefficient of lift due to pitch rate 
coefficient of lift due to speed brake deflection 
coefficient of rolling moment due to roll rate 
coefficient of rolling moment due to yaw rate 
total coefficient of pitching moment 
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I 

I 

CMA coefficient of pitching moment due to angle of attack 

CMAD coefficient of pitching moment due to angle-of-attack rate 

CMDE coefficient of pitching moment due to symmetric elevator deflection 

CMO pitching moment coefficient at zero angle of attack 

CMQ coefficient of pitching moment due to pitch rate 

CMSB coefficient of pitching moment due to speed brake deflection 

CN total coefficient of yawing moment 

i 

1 CNB coefficient of yawing moment due to sideslip 

CNDA coefficient of yawing moment due to aileron deflection 

CNDR coefficient of yawing moment due to rudder deflection 

CNDT coefficient of yawing moment due to differential elevator deflection 

CNP coefficient of yawing moment due to roll rate 

CNR coefficient of yawing moment due to yaw rate 

CY total coefficient of sideforce 

CYB coefficient of sideforce due to sideslip 

CYDA coefficient of sideforce due to aileron deflection 

CYDR coefficient of sideforce due to rudder deflection 

CYDT coefficient of sideforce due to differential elevator deflection 

DAS lateral trim parameter 

DC surface deflection and thrust control array 

DELX displacement of the aerodynamic reference point along the x body axis from 

the vehicle center of gravity 

DELY displacement of the aerodynamic reference point along the y body axis from 

the vehicle center of gravity 

DELZ displacement of the aerodynamic reference point along the z body axis from 

the vehicle center of gravity 

DES longitudinal trim parameter 
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DRS 

DXTHRS 

EIX 

ENGOMG 

GMA 

H 

HDOT 

NUMSAT 

NUMSRF 

P 

PDOT 

PHI 

PHIDOT 

PSI 

PSIDOT 

Q 

QBAR 

QDOT 

R 

RDOT 

S 

T 

TDOT 

THA 

THADOT 

THRSTX 

THRUST 


directional trim parameter 

distance between the center of gravity of the engine and the 
thrust point 

rotational inertia of the engine 
rotational velocity of the engine 
flightpath angle 
altitude 

time rate of change of altitude 
number of states in the state vector 
number of recognizable control names 
roll rate 

time rate of change of roll rate 
roll angle 

time rate of change of roll angle 
heading angle 

time rate of change of heading angle 
pitch rate 
dynamic pressure 

time rate of change of pitch rate 
yaw rate 

time rate of change of yaw rate 

wing area 

time 

time rate of change of time 
pitch angle 

time rate of change of pitch angle 

thrust trim parameter 

thrust generated by each engine 
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> 


TLOCAT 


TVANXY 

TVANXZ 

UB 

V 

VB 

VCAS 

VDOT 

VEAS 

WB 

X 


XDOT 

XYANGL 

XZANGL 

Y 

YDOT 


location of the engine in the x, y, and z axes from the center of gravity 
orientation of the thrust vector in the x-y engine axis plane 
orientation of the thrust vector in the x-z engine axis plane 
velocity along the x body axis 
velocity 

velocity along the y body axis 
calibrated airspeed 

time rate of change of total vehicle velocity 

equivalent airspeed 

velocity along the z body axis 

position north from an arbitrary reference point 
time rate of change of north-south position 
orientation of thrust vector in x-y body axis plane 
orientation of thrust vector in x-z body axis plane 
position east from an arbitrary reference point 
time rate of change of east-west position 


PROGRAM OVERVIEW 

The program LINEAR numerically determines a linear system model using nonlinear 
equations of motion and a user-supplied nonlinear aerodynamic model. LINEAR is also 
capable of extracting linearized gross engine effects, such as net thrust, torque, 
and gyroscopic effects and including these effects in the linear system model. The 
point at which this linear system model is defined is determined either by specifying 
the state and control variables or by selecting an analysis point on a trajectory, 
selecting a trim option, and allowing the program to determine the control variables 
and remaining state variables to satisfy the trim option selected. 

Because the program is designed to satisfy the needs of a broad class of users, 
a wide variety of options has been provided. Perhaps the most important of these 
options are those that allow user specification of the state, control, and obser- 
vation variables to be included in the linear model derived by LINEAR. 

Within the program, the nonlinear equations of motion include 12 states repre- 
senting a rigid aircraft flying in a stationary atmosphere over a flat, nonrotating 
earth. Thus, internally, the state vector x is computed as 

x = [p q r V a £ <j> 0 i|;hxyJ T 
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The nonlinear equations used to determine the derivatives of the states are pre- 
sented in the following section. The internal control vector u can contain up to 
30 controls. The internal observation vector y contains 120 variables including the 
state variables, the time derivatives of the state variables, the control variables, 
and a variety of other parameters of interest. Thus, within the program. 



where 


Y 1 


L a x 

a y a z a. 

nx a ny 

a a a • a • a • a • 

nz n nx,i ny,i nz,i n,i 

Y2 

= 

[ a 

Re Re ' M 

q q c 

Pa qc/Pa Pt T T t v e V C ] T 

Y3 

= 

[y 

• •• 
fpa y h 

h/57.3 

f 



r 

„ iT 



Y4 


[e s 

p s J 



Y5 

= 

[ L 

D N a] T 





r 

• • 

• iT 


Y6 


l u 

V W U V 

wj 


Y7 

= 

[ a ,i 

*.i h ,i 

■Uf 



and 

Y8 = [t p s q s r sJ T 


The equations defining these quantities are presented in the Observation Equations 
section. 


From the internal formulation of the 
state, control, and observation variables, 
the user must select the specific vari- 
ables desired in the output linear model 
(described in the Selection of State, Con- 
trol, and Observation Variables section). 
Internally, the program uses a 12-state 
vector, 120-variable observation vector, 
and a 30 -parameter control vector. These 
variables can be selected to specify the 
formulation most suited for the specific 
application. The order and number of 
parameters in the output model are com- 
pletely under user control. Figure 1 
shows the selection of the variables in 
the state vector for a requested linear 
model. From the internal formulation on 


Output model Internal 

parameters parameters 



P 

q 

r 

V 

a 

P 

<{> 

e 

v 

h 

x 

y 


Figure 1. Selection of state vari- 
ables for linear model. 
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the right, the requested model is constructed, and the linear system matrices are 
selected in accordance with the user specification of the state, control, and obser- 
vation variables. 

The linear model derived by LINEAR is determined at an analysis point. LINEAR 
allows this analysis point to be defined as a true steady-state condition on a spec- 
ified trajectory (a point at which the rotational and translational accelerations 
are zero) or a totally arbitrary state on a trajectory. The program LINEAR provides 
the user with several options described in detail in the Analysis Point Definition 
section of this report. These analysis point definition options allow the user to 
trim the aircraft in wings-level flight, pushovers, pullups, level turns, or zero- 
sideslip maneuvers. Also included is a nontrimming option in which the user defines 
a totally arbitrary condition about which the linear model is to be derived. 

The linear system matrices are determined by numerical perturbation and are the 
first-order terms of a Taylor series expansion about the analysis point as described 
in the Linear Models section. The formulation of the output system model is under 
user control. The user can select state equation matrices corresponding to either 
the standard formulation of the state equation 

• 

X = Ax + Bu 


or the generalized equation 


Cx = A'x + B'u 


The observation matrices can be selected from either of two formulations 
corresponding to the standard equation 

y = Hx + Fu 


or the generalized equation 


y=H'x+Gx+F , u 

In addition to the linear system matrices, LINEAR also computes the nondimen- 
sional stability and control derivatives at the analysis point. These derivatives 
are discussed in the Aerodynamic Model section. 

Aircraft geometry and mass properties and analysis point definition are input 
interactively to the program or are read from ASCII files of 80-character records 
specified by the user. The state, control, and observation vectors are also de- 
fined in either manner. The details of data input are discussed in the Data In- 
put section. 

The output of LINEAR consists of the following four files: one is intended to 

be used with follow-on design and analysis programs, two document the options and 
analysis points selected by the user, and one contains stability and control deriva- 
tives extracted at the analysis point. The first of these files contains the state, 
control, and observation vector name lists and values, and the state space system 
matrices. The second file contains all of the trim information, the control, state. 
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and observation values, the trim control input positions, and the nondimensional 
stability and control derivative data. It will also contain the output of the 
equations of motion if the program fails to attain a trimmed condition. The third 
file contains all of the information in the second file plus the state space system 
model matrices. The fourth file contains the stability and control derivatives 
extracted at the analysis point in a form suitable to be read back into LINEAR as 
a linear aerodynamic model. 

To execute LINEAR, five user-supplied subroutines are required. These routines, 
discussed in the User-Supplied Subroutines section, define the nonlinear aerodynamic 
model, the gross engine model, the gearing between the LINEAR trim inputs and the 
surfaces modeled in the aerodynamic model, and a model of the mass and geometry 
properties of the aircraft. The control model, shown in figure 2, defines how the 
LINEAR trim inputs will be connected to the surface models and allows schedules and 
nonstandard trimming schemes to be used. The last feature is particularly important 
for asymmetrical aircraft. 


Inputs from 
/CTPARM/ 


DES 
DAS- 
DRS ■ 
THRSTX • 


UCNTRL 

(Gearing) 


Outputs to 
/CONTRL/ 

DC(1) 

*- DC(2) 


DC ( 30 ) 


(Pilot stick, pedal, 
and throttle) 


(Surface deflections 
and power level setting) 


Figure 2. Inputs and outputs to the 
user-supplied subroutine UCNTRL . 


EQUATIONS OF MOTION 


The nonlinear equations of motion used in the linearization program are general 
six-degree-of-freedom equations representing the flight dynamics of a rigid aircraft 
flying in a stationary atmosphere over a flat nonrotating earth. The assumption of 
nonzero forward motion also is included in these equations. Because of this assump- 
tion, these equations are invalid for vertical takeoff and landing or hover. These 
equations contain no assumptions of either symmetric mass distribution or aerody- 
namic properties and are applicable to asymmetric aircraft (such as oblique-wing 
aircraft) as well as to conventional symmetric aircraft. These equations of motion 
were derived by Etkin (1972), and the derivation is detailed in the NASA Reference 
Publication, "Derivation and Definition of a Linear Aircraft Model," by Eugene L. 
Duke, Robert F. Antoniewicz, and Keith D. Krambeer. 

The following equations for rotational acceleration are used for analysis point 
definition: 
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where 


p = [ ( EL ) 1 1 + (EM)l2 + (E1OI3 “ P^dxz I 2 “ ^xy 3 ^ 

+ PqUxz 1 ! “ ^yz *2 “ ^z^ 3 ^ ~ P r ^xy ^1 + *^y ^2 ~ ^yz^ 3 ^ 
+ q 2 (Iyz I 1 ” I xy I 3 ^ ~ qr(D x I 1 - I xy l 2 + I xz I 3 ) 

- r 2 (I yz Ii - I xz I 2 )]/det I 

q = [( EL)I 2 + (XM)I 4 + (ZN)I 5 - P 2 (I xz I 4 - I X y I 5 ) 

+ Pq( I xz I 2 “ -*-yz^4 ” *^z^5^ ~ P r ^ I xy^2 + °y I 4 ^yz^S ^ 

+ q 2 (Iyz I 2 “ I xy I 5^ “ < J r ^ D x I 2 “ I xy I 4 + ^xz^^ 

- r 2 (I yz I 2 - I xz I 4 )]/det I 

r = [(EDI3 + (EM)l5 + (EN)Ig - P 2 (I xz l5 “ I xy I 6^ 

+ pq(I xz l 3 - I yz l 5 - D z Ig) - pr(I xy l 3 + D y l 5 “ lyz^ 6 ^ 

+ q 2 (Iyz I 3 “ I xy I 6^ ” < 5 r ^ D x I 3 - I xy I 5 + 

- r 2 (I yz l3 - I xz l5 ) ] /det I 


det I 

— 

■^■x^y-'-z " 2 ^xy^ 

II 

= 

lylz - ^- 2 yz 

12 

= 

^xy I z + ^yz^xz 

*3 

= 

^■xy^yz + ^y^xz 

14 

= 

•^x^z “ l2 xz 

15 

= 

-^x-^yz + ^xy-^xz 

16 

= 

I x I y ” l2 xy 

Dx 

= 

Iz ~ ^y 

>1 

Q 

= 

*x - I Z 

D z 

= 

ly " *x 


I x l2 yz “ ^y-^'^xz 


I z l2 xy 


Here, the body axis rates are designated p, q, and r, corresponding to roll rate, 
pitch rate, and yaw rate. The total moments about the x, y, and z body axes are 
designated EL, EM, and EN, respectively. These total moments are the sums of all 
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aerodynamic moments and powerplant- induced moments due to thrust asymmetries and 
gyroscopic torques. The equations used to determine the change in moment coef- 
ficients due to the noncoincidence of the vehicle center of gravity and the refer- 
ence point of the nonlinear aerodynamic model are derived in appendix A. The equa 
tions defining the engine torque and gyroscopic contributions to the total moments 
are derived in appendix B. The body axis moments of inertia are designated I x , l y 
and I z . The products of inertia are designated I xy , I xz , and I yz . These moments 
and products of inertia are elements of the inertia tensor I, defined as 


Ix 

“Ixy 

-Ixz 

“Ixy 

X Y 

“lyz 

“Ixz 

N 

>i 

H 

1 

Iz 


To derive the state equation matrices for the generalized formulation 

Cx = Ax + Bu 

the rotational accelerations are cast in a de coupled-axes formulation. The equa- 
tions used to derive the linearized matrices are 


1 Ixy "Ixz 


p 

Ix Ix 



“Ixy “lyz 

■1 


• 

H 

H 


q 

“Ixz -lyz 


• 

r 

Iz Iz 


L. J 


EM Ix Ixy lyz 

rp — + rq L - pq -f— 

iy J-y -Ly J-y 


+ ( r* 


2 , X xz Iz 

p**) J + pr — 

y y 


EN Ix Ixz lyz , . o o . Ixy Iy 

— + qp •= qr + pr -f— + (p^ - q^ ) — -i - pq 

A z A z A z L z L z L z 


The translational acceleration equations used in the program LINEAR for both 
analysis point definition and perturbation are 

V = [— D cos B + Y sin 0 + X T cos a cos 0 + Y T sin 0 + Z T sin a cos 0 

- mg(sin G cos a cos 0 - cos 0 sin <j> sin 0 

- cos 0 cos d> sin a cos 0)]/m 
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a = [-L + Zip cos a - X T sin a + mg(cos 6 cos <j> cos a 
+ sin 0 sin a)]/Vm cos 0 

+ q - tan 0 (p cos a + r sin a) 

0 = [D sin 0 + Y cos 0 - X T cos a sin 0 + Y T cos 0 - Z T sin a sin 0 

+ mg(sin 0 cos a sin 0 + cos 0 sin <p cos 0 

- cos 0 cos <(> sin a sin 0)]/Vm 
+ p sin a - r cos a 

The equations used to define the vehicle attitude rates are 

= p + q sin <j> tan 0 + r cos <p tan 0 

0 = q cos <j> - r sin (p 

ip = q sin <p sec 0 + r cos <p sec 0 

The equations used to define the earth-relative velocities are 

h = v(cos 0 cos a sin 0 - sin 0 sin <p cos 0 - cos 0 sin a cos <j> cos 0) 

x = V[cos 0 cos a cos 0 cos ip + sin 0 (sin <p sin 0 cos <p - cos <p sin ip) 

+ cos 0 sin a (cos <f> sin 0 cos tp + sin <j> sin ip)] 

y = V[cos 0 cos a cos 0 sin tp + sin 0 (sin <p sin 0 sin ip + cos <p cos ip) 

+ cos 0 sin a (cos <p sin 0 sin ip - sin <p cos \p)] 

OBSERVATION EQUATIONS 

The user-selectable observation variables computed in LINEAR represent a broad 
class of parameters useful for vehicle analysis and control design problems. These 
variables include the state, time derivatives of state, and control variables. Also 
included are air data parameters, accelerations, flightpath terms, and other miscel- 
laneous parameters. The equations used to calculate those parameters are derived 
from a number of sources (Clancy, 1975; Dommasch and others, 1967; Etkin, 1972; 
Gainer and Hoffman, 1972; Gracey, 1980). Implicit in many of these observation 
equations is an atmospheric model. The model included in LINEAR is derived from 
the U.S. Standard Atmosphere (1962). 

The vehicle body axis accelerations constitute the set of observation variables 
that, except for state variables themselves, are most commonly used in the aircraft 
control analysis and design problem. These accelerations are measured in g units 
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and are derived directly from the body axis forces defined in the previous section 
for translational acceleration. The body axis acceleration equations used in LINEAR 
are 



The equations 
are 


(Xip - D cos a + L sin a - gm sin 0)/g Q m 

(Y T + Y + gm cos 0 sin <|>)/g o m 

(Zfp - D sin a - L cos a + gm cos 0 cos <j>)/g Q m 

for body axis accelerometers that are at vehicle center of gravity 


a n ^ = (Xip - D cos a + L sin a)/g Q m 


a n y = (* T + Y >/9 0 m 


a n ^ = ( Zip - D sin a - L cos a)/g Q m 


*n 


(-Z T + D sin a + L cos a)/g m 


For orthogonal accelerometers that are aligned with the vehicle body axes but are 
not at vehicle center of gravity, the following equations apply: 

a n x ,i = a n x " [( <1 2 + r2 >*x " <PS " *>Yx " ( P r + < 3 )z x^/g 0 
a n y/i = a n y + [(Pq + r ) x y ~ (P 2 + r 2 )y y + (qr - p)z y ]/g o 
a n z , i = a n z + t ( P r ” <3) x z + (< 3 r + P^z “ <<1 2 + P 2 )z z J/g o 
a n,i = a n “ [(pr - q)x z + (qr + p)y z - (q2 + p 2 )z z ]/g Q 


where the subscripts x, y, and z refer to the x, y, and z body axes, respectively, 
and the symbols x, y, and z refer to the x, y, and z body axis locations of the sen- 
sors relative to the vehicle center of gravity. The symbol g Q is the acceleration 

due to gravity at sea level. Also included in the set of acceleration parameters 
is the load factor n = L/W, where L is the total aerodynamic lift and W is the 
vehicle weight. 


The air data parameters having the greatest application to aircraft dynamics and 
control problems are the sensed parameters and the reference and scaling parameters. 
The sensed parameters are impact pressure q c , ambient or free-stream pressure p a , 
total pressure p t , ambient or free-stream temperature T, and total temperature T t . 

The selected reference and scaling parameters are Mach number M, dynamic pressure q, 
speed of sound a, Reynolds number Re, Reynolds number per unit length Re', and the 
Mach meter calibration ratio q c /p a . These quantities are defined as 


a 


( 


1 .4 


Po 

PoT 0 


T ) 


1/2 
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M 


V 

a 


R c 


= P v£ 


R ; --f 


pV 2 

q = ■— 



- 1.0 


(1.0+ 0.2M2) 3 * 5 - 1 .( 


f 

[i.2m 2(_ \ 2 ’ 5 ] . 1e0 

k |_ \ 5.6M 2 - 0.8/ J 


P t = *a + 


T t = T( 1.0 + 0.2M 2 ) 


(M < 1.0) 
(M > 1.0) 

(M < 1.0) 
(M > 1.0) 


where p is the density of the air, p is the coefficient of viscosity, and the sub- 
script o refers to sea level, standard day conditions. Free-stream pressure, free- 
streara temperature, and the coefficient of viscosity are derived from the U.S. Stand- 
ard Atmosphere (1962). 

Also included in the air data calculations are two velocities: equivalent air- 

speed V e and calibrated airspeed V c , both computed in knots. The calculations assume 
that internal units are in the English system. The equation used for equivalent air- 
speed is 


v e = 17. 17 Vq 


which is derived from the definition of equivalent airspeed. 



where p Q = 0.002378 slug/ft 2 and V e is converted from feet per second to knots. Cal- 
ibrated airspeed is derived from the following definition of impact pressure: 
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[( Po 2 \ 3 * 5 

mP o[{'' 0 + K*o W °) ' ’J 


(V c < a Q ) 


= 1 - 2 r- p 


v c x2 


5.76 


2.5 


5.6 - 0.8 (a 0 /V c ) J 
For the case where V c < a Q , the equation for V c is 


" P Q ( v c > a o> 


V c = 1479.116 


qc 


2/7 

+ 1.0) - 1.0 (V c < a D ) 


Calibrated airspeed is found using an iterative process for the case where V c > a D . 
The equation 


V c = 582.95174 [( — + 1 . o\ fl.O — (V c > a 0 } 

VV P o / |_ 7.0 (V c /a 0 ) 2 j 

is executed until the change in V c from one iteration to the next is less than 
0.001 knots. 


Also included in the observation variables are the f lightpath-related parameters 
(described in app. C) , including flightpath angle y, flightpath acceleration fpa, 

'* • 

vertical acceleration h, flightpath angle rate y, and (for lack of a better category 

in which to place it) scaled altitude rate h/57.3. The equations used to determine 
these quantities are 


y = sm 


-ii) 


£ V 

fpa = — 

9o 


h = a x sin 0 - a y sin 0 cos 0 - a z cos <)> cos 0 


y = 


V h - hV 

v Vv 2 - h 2 


Two energy-related terms are included in the observation variables: specific 

energy E s , and specific power P s , defined as 


E 


s 


h 


+ 


V 2 

2g 


p 


s ~ 


h 


+ 


VV 
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The set of observation variables available in LINEAR also includes four force 
parameters: total aerodynamic lift L, total aerodynamic drag D, total aerodynamic 

normal force N, and total aerodynamic axial force A. These quantities are defined as 


L = qSC L 


D — qSCp 

N = L cos a + D sin a 
A = -L sin a + D cos a 

where Cp and Cl are coefficients of drag and lift, respectively. 

Six body axis rates and accelerations are available as observation variables. 
These include the x body axis rate u, the y body axis rate v, and the z body axis 

• • • 

rate w. Also included are the time derivatives of these quantities, u, v, and w. 
The equations defining these quantities are 

u = V cos a cos 0 

v = V sin $ 

w = V sin a cos (3 


- gm sin 0 - D cos 


a + L sin 


+ rv - qw 


i-(is 

• / y T + 9 111 cos ® sin <t> + Y\ 

v \ m / 

• ^ Z T + 9 111 cos 8 cos <j> ~ D sin a - L cos 


+ pw - ru 


+ qu - pv 


The final set of observation variables available in LINEAR is a miscellaneous 
collection of other parameters of interest in analysis and design problems. The 
first group consists of measurements from sensors not located at the vehicle center 
of gravity. These represent angle of attack a f i, angle of sideslip 0 altitude 

h f i, and altitude rate h^ measurements displaced from the center of gravity by some 
x, y, and z body axis distances. The equations used to compute these quantities are 

3 * T .JSSC 

V 


«.i - ° - ( 2L ^) 


h f i = h + x sin 9 - y sin <£ cos 8 - z cos <J> cos 0 
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h^ = h + 0(x cos 0 + y sin <f> sin 0 + z cos <J> sin 0) 


- 4> (y cos <j) cos 0 - z sin $ cos 0) 

The remaining miscellaneous parameters are total angular momentum T, stability axis 
roll rate p s , stability axis pitch rate q s , and stability axis yaw rate r s , defined 
as 


T = J^xP 2 - 2I xy pq - 2I xz pr + I y q2 - 2I yz qr + I z r2 ) 

p s = p cos a + r sin a 
^s = 9 

r s = -p sin a + r cos a 

SELECTION OF STATE, CONTROL, AND OBSERVATION VARIABLES 


The equations in the two preceding sections define the state and observation 
variables available in LINEAR. The control variables are defined by the user in the 
input file. Internally, the program uses a 12-state model, a 120-variable obser- 
vation vector, and a 30-parameter control vector. These variables can be selected 
to specify the formulation most suited for the specific application. The order and 
number of parameters in the output model are completely under user control. Figure 1 
shows the selection of variables for the state vector of the output model. However, 
it should be noted that no model order reduction is attempted. Elements of the 
matrices in the internal formulation are simply selected and reordered in the for- 
mulation specified by the user. 

Specific state, control, and observation variables for the formulation of the 
output matrices are selected by alphanumeric descriptors in the input file. The use 
of these alphanumeric descriptors is described in the Data Input section. Appendix C 
lists the observation variables and their alphanumeric descriptors. Appendix D lists 
the state variables and their alphanumeric descriptors. The alphanumeric descriptors 
for the selection of control parameters to be included in the observation vector are 
the control variable names defined by the user in the input file, as described in the 
Input Files section. 


LINEAR MODELS 


The linearized system matrices computed by LINEAR are the first-order terms of 
a Taylor series expansion about the analysis point (Dieudonne, 1978; Kwakernaak and 
Sivan, 1972; NASA RP by Duke, Antoniewicz, and Krambeer) and are assumed to result in 
a time-invariant linear system. The validity of this assumption is discussed in the 
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Analysis Point Definition section. The technique used to obtain these matrices 
numerically is a simple approximation to the partial derivative, that is, 

3f f ( x 0 + Ax) - f(x Q - Ax) 

3 X 2 Ax 

where f is a general function of x, an arbitrary independent variable. The Ax may be 
set by the user, but defaults to 0.001 for all state and control parameters, with the 
single exception of velocity V, where Ax is multiplied by a, the speed of sound, to 
obtain a reasonable perturbation size. 

From the generalized nonlinear state equation 

Tx = f(x, x, u) 


and the observation equation 


y = g(x, x, u) 

the program determines the linearized matrices for the generalized formulation of the 
system: 


C <5x = A' 6x + B ' 6u 

6y = H 1 6x + G 6x + F* 6u 


where 


C 


T 


3f 

3x 


A' 


3f 

3x 


B' 


3f 

3u 


K' = |2 

3x 


G 


ia 

3x 


F' 


ia 

3u 


23 



with all derivatives evaluated along the nominal trajectory defined by the analysis 

point (x D , x Q/ u 0 ) . The state, time derivative of state, and control vectors can be 
expressed as small perturbations about the nominal trajectory, so that 

x = x Q + 6x 

• • • 

x = x 0 + 6x 

u = u 0 + 6u 

In addition to the matrices for this generalized system, the user has the option 
of requesting linearized matrices for a standard formulation of the system: 


where 


6x = A 6x + B 6u 
6y = H 6x + F Su 


A 


T 



3f 

3x 


B 


T 


_3f 

3x 


-1 


3f 

3u 


H 



3g 

3x 


T 


3x 


-1 


3f 

3x 



3g 

3x 


T 


_3f 

3x 


-1 


3f 

3u 


with all derivatives evaluated along the nominal trajectory defined by the analysis 
point (x Q , x 0 , u Q ) . 

LINEAR also provides two nonstandard matrices, D and E, in the equations 

x = Ax + Bu + Dv 


y = Hx + Fu + Ev 


or D' and E' in the equations 


Cx = A ' x + B'u + D’v 
y = H'x + Gx + F'u + E'v 
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These dynamic interaction matrices include the effect of external forces and moments 
acting on the vehicle. The components of the dynamic interaction vector v are incre- 
mental body axis forces and moments: 


'fix’ 

6Y 



6M 

L<5n_ 


Thus, 



E' 


= is 

9v 


and 



9f "I -1 _9f 
9xJ 9v 


E 


= is + is 

3v 9x 





9f 

9v 


These matrices allow the effects of unusual subsystems or control effectors to be 
easily included in the vehicle dynamics. 


The default output matrices for LINEAR are those for the standard system for- 
mulation. However, the user can select matrices for either generalized or standard 
state and observation equations in any combination. Internally, the matrices are 
computed for the generalized system formulation and then combined appropriately to 
accommodate the system formulation requested by the user. 


ANALYSIS POINT DEFINITION 


The point at which the nonlinear system equations are linearized is referred to 
as the analysis point. This can represent a true steady-state condition on the spec- 
ified trajectory (a point at which the rotational and translational accelerations are 
zero? Perkins and Hage, 1949; Thelander, 1965) or a totally arbitrary state on a tra- 
jectory. LINEAR allows the user to select from a variety of analysis points. These 
analysis points are referred to as trim conditions, and several options are available 
to the user. The arbitrary state and control option is designated NOTRIM, and in 
selecting this option the user must specify all nonzero state and control variables. 
For the equilibrium conditions, the user specifies a minimum number of parameters, 
and the program numerically determines required state and control variables to force 
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the rotational and translational accelerations to zero. The analysis point options 
are described in detail in the following subsections. 

For all the analysis point definition options, any state or control parameter 
may be input by the user. Those state variables not required to define the analysis 
point are used as initial estimates for the calculation of the state and control 
conditions that result in zero rotational and translational accelerations. As each 
state variable is read into LINEAR, the name is compared to the list of alternative 
state variable names listed in appendix D. All state variables except velocity 
must be specified according to this list. Velocity can also be defined by specify- 
ing Mach number (see alternative observation variable names in appendix C). Appen- 
dix E lists analysis point definition identifiers that are recognized by LINEAR. 

It should be noted that the option of allowing the user to linearize the system 
equations about a nonequilibrium condition raises theoretical issues (beyond the 
scope of this report) the potential user should be aware of. Although all the anal- 
ysis point definition options provided in LINEAR have been found to be useful in 
the analysis of vehicle dynamics, not all the linear models derived about these anal- 
ysis points result in the time-invariant systems assumed in this report. However, 
the results of the linearization provided by LINEAR do give the appearance of being 
time invariant. 

The linearization process as defined in this report is always valid for some time 
interval beyond the point in the trajectory about which the linearization is done. 
However, for the resultant system to be truly time invariant, the vehicle must be in 
a sustainable, steady-state flight condition. This requirement is something more 

than merely a trim requirement, which is typically represented as x(t) = 0, indi- 
cating that for trim, all the time derivatives of the state variables must be zero. 
(This is not the case, however: Trim is achieved when the acceleration-like terms 

are identically zero; no constraints need to be placed on the velocity-like terms in 

• •••••• 

x. Thus, for the model used in LINEAR, only p, q, r, V, a, and 0 must be zero to 
satisfy the trim condition. ) The trim condition is achieved for the straight-and- 
level, pushover-pul lup, level turn, thrust-stabilized turn, and beta trim options 
described in the following sections. In general, the untrimmed and specific power 
analysis point definition options do not result in a trim condition. 

Of these analysis point options resulting in a trim condition, only the straight- 
and-level and level turn options force the model to represent sustainable flight con- 
ditions. In fact, only in the special case where the flightpath angle is zero does a 
truly sustainable trim occur. 

As previously stated, the linearization of a nonlinear model and its represen- 
tation as a time-invariant system are always valid for some time interval beyond the 
analysis point on the trajectory. This time interval is determined by several fac- 
tors (such as trim and sustainable flight conditions) and ultimately by accuracy 
requirements placed on the representation. Thus, in using the linear models provided 
by this program, the user should exercise some caution. 


Untrimmed 

For the untrimmed analysis point, the user specifies all state and control var- 
iables that are to be set at some value other than zero. The number of state var- 


26 



iables specified is entirely at the user's discretion. If any of the control var- 
iables are to be nonzero, the user must specify the control parameter and its value. 
The untrimmed option allows the user to analyze the vehicle dynamics at any flight 
condition, including transitory conditions. 


Straight-and-Level Trim 

The straight-and-level analysis points available in LINEAR are in fact wings- 
level, constant-f lightpath-angle trims. Both options available for straight-and- 
level trim allow the user to specify either a flightpath angle or an altitude rate. 
However, since the default value for these terms is zero, the default for both types 
of straight-and-level trim is wings-level, horizontal flight. 

The two options available for straight-and-level trim require the user to specify 
altitude and either an angle of attack or a Mach number. If a specific angle of 
attack and altitude combination is desired, the user selects the "Mach-trim" option, 
which determines the velocity required for the requested trajectory. Likewise, the 
"alpha-trim" option allows the user to specify Mach number and altitude, and the trim 
routines in LINEAR determine the angle of attack needed for the requested trajectory. 

The trim condition for both straight-and-level options is determined within the 
following constraints: 


p = q = r = 0 
<p - 0 

The trim surface positions, thrust, angle of sideslip, and either velocity or angle 
of attack are determined by numerically solving the nonlinear equations for the 
translational and rotational acceleration. Pitch attitude 0 is determined by itera- 
tive solution of the altitude rate equation. 


Pushover-Pullup 

The pushover-pul lup analysis points result in wings-level flight at n * 1 . For 

n > 1, the analysis point is the minimum altitude point of a pullup when h = 0. For 

n < 1, this trim results in a pushover at the maximum altitude with h = 0. There 
are two options available for the pushover-pullup analysis point definition: (1) the 

"alpha-trim" option in which angle of attack is determined from the specified alti- 
tude, Mach number, and load factor, and (2) the "load-factor-trim" option in which 
angle of attack, altitude, and Mach number are specified, and load factor is deter- 
mined according to the constraint equations. 

The analysis point is determined at the specified conditions, subject to the fol- 
lowing contraints: 

p = r = 0 

<1 = r, n [mg(n - cos (8 - a)) - Z T cos a + Xm sin a] 

Vm cos p A x 

<p = 0 
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The expression for q is derived from the a equation by setting a = 0 and <j> = 0; 0 is 

• 

derived from the h equation. The trim surface positions, thrust, angle of sideslip, 
and either angle of attack or load factor are determined by numerically solving the 
nonlinear equations for the translational and rotational accelerations. 


Level Turn 

The level turn analysis point definition options result in nonwings-level, 
constant-turn-rate flight at load factors greater than one. The vehicle model is 
assumed to have sufficient excess thrust to trim at the condition specified. If 
thrust is not sufficient, trim will not result, and the analysis point thus defined 
will have a nonzero (in fact, negative) velocity rate. 

The level trim options available in LINEAR require the specification of an alti- 
tude and a Mach number. The user can then use either angle of attack or load factor 
to define the desired flight condition. These two options are referred to as "alpha- 
trim" and "load-factor-trim," respectively. For either option, the user may also 
request a specific flightpath angle or altitude rate. Thus, these analysis point 
definitions may result in ascending or descending spirals, although the default is 
for the constant-altitude turn. 

The constraint equations for the coordinated level turn analysis point defini- 
tions are derived by Chen (1981), and Chen and Jeske (1981). Using the requested 
load factor, the tilt angle of acceleration normal to the flightpath from the ver- 
tical plane <{>l is calculated using the equation 

= ±tan-1 [ (n 2 - cos 2 y)„ 1/2 1 
^ |_ cos y J 

where the positive sign is used for a right turn and the negative sign is used for a 
left turn. From <f> L , turn rate can be calculated as 

^ = V tan ^ 


Using these two definitions, the state variables can be determined: 
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0 = sin“1(-p/ty) 


<(> = tan -1 (q/r) 

The trim surface positions, thrust, angle of sideslip, and either angle of attack or 
load factor are determined by numerically solving the nonlinear equations for the 
translational and rotational accelerations. 


Thrust-Stabilized Turn 

The thrust-stabilized turn analysis point definition results in a constant throt- 
tle nonwings-level turn with a nonzero altitude rate. The two options available in 
LINEAR are "alpha-trim" and "load-factor-trim." These options allow the user to spe- 
cify either the angle of attack or the load factor for the analysis point. The alti- 
tude and Mach number at the analysis point must be specified for both options. The 
user also must specify the value of the thrust trim parameter by assigning a value to 
the variable THRSTX in the input file after the trim has been selected. 

The constraint equations for the thrust-stabilized turn are the same as those for 
the level turn. However, for this analysis point definition, flightpath angle is 
determined by LINEAR. 


Beta Trim 

The beta trim analysis point definition results in nonwings-level, horizontal 

flight with the heading rate \p equal to zero at a user-specified Mach number, alti- 
tude, and angle of sideslip. This trim option is nominally at 1 g, but as 3 and <j> 
vary from zero, normal acceleration decreases and lateral acceleration increases. 

For an aerodynamically symmetric aircraft, a trim to 3 = 0 using the beta trim option 
results in the same trimmed condition as the straight-and- level trim. However, for 
an aerodynamically asymmetric aircraft, such as an oblique-wing vehicle, the two trim 
options are not equivalent. 

The constraint equations used with the beta trim option are 

p = q = r = 0 

The trim surface positions, thrust, angle of attack, and bank angle are determined by 
numerically solving the nonlinear equations for translational and rotational acceler- 
ation. Pitch attitude 0 is derived from the equation for flightpath angle y with 
Y = 0 and is defined as 

q _ t an -1 / sin 3 sin + cos 3 sin g \ 

\ cos 3 cos a ) 


Specific Power 

The specific power analysis point definition results in a level turn at a user- 
specified Mach number, altitude, thrust trim parameter, and specific power. Unlike 
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the other trim options provided in LINEAR, the specific power option does not, in 

general, attempt to achieve the zero velocity rate V. In fact, because the altitude 

* 

rate h = 0 and specific power are defined by 


P 


s 



the resultant velocity rate will be 


V 


1*1 

V 


However, the other acceleration-like terms (p, q, r, a, and |3) will be zero if the 
requested analysis point is achieved. 

The constraint equations used with the specific power analysis point definition 
option can be derived from the load factor tilt angle equation used in the level turn 
analysis point definition with y = 0: 

<j)L = ±tan -1 (n 2 - 1) 1//2 

(where the positive sign is used for a right turn and the negative sign is used for a 
left turn). 


^ ~ y tan <Pl 
q = sin <f>L cos g 
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s 


2 

tan <{> L cos g 


p s = -q tan g 
p = p s cos u - r s sin a 
r = p s sin a + r s cos a 
0 = sin“l / - 

\ V 

<t> = tan -- ' SL 
T r 

The analysis point surface positions, load factor, angle of attack, and angle of 
sideslip are determined by numerically solving the nonlinear equations for transla- 
tional and rotational acceleration. 


30 


NONDIMENSIONAL STABILITY AND CONTROL DERIVATIVES 


The nondimensional stability and control derivatives computed by LINEAR from the 
nonlinear aerodynamic model assume broadly formulated linear aerodynamic equations. 
These equations include effects of what are normally considered exclusively lateral- 
directional parameters in the longitudinal force and moment coefficient equations, 
and, conversely, effects of longitudinal parameters in the lateral-directional equa- 
tions. The reason for this is two-fold: application to a larger class of vehicle 

types, such as oblique-wing aircraft, and computational ease. 

The nondimensional stability and control derivatives assume the following equa- 
tions for the aerodynamic force and moment coefficients: 

n 

= c £o + c A a a + C *0 0 + C * h 6h + C *V 6v + 2 c *6- 6 i 

i=1 1 


+ c 




II 

Cm - CmQ + Cjj^a + C^B + C m ^ 6h + C,^ 6V + ^ C m< ^ ^i 


+ ^pP + C m q q + C m r r + + ^B^ 
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n a a + c ng3 + C nh 6h + C nv 6V + J, C n6 . 6 ± 


+ c n p P + C n q + C nr r + C n *u + c *0 


C L = c Lo + c L a “ + C L g3 + C Lh 5h + C Lv 6v + £ C L 6i 

i=1 1 


+ C Lp p + C Lq q + C Lr r + C L -a + C l *B 


C D = C Do + C Da a + C D3 B + C Dh 5h + C Dv 6V + X C D6 . 6 A 
+ C Dp p + C Dq q + C Dr f + C D *a + C D *B 
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+ C Y p P + C Yq-5 + C Y r r + C Ya a + C Y 3 ^ 
where the stability and control derivatives have the usual meaning, 

3C e 

c e = t — 
b x 9x 


with C e being an arbitrary force or moment coefficient and x an arbitrary nondimen- 
sional variable. The rotational terms in the equations are nondimensional versions 
of the corresponding state variable with 
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The 6^ in the summations are the n control variables defined by the user. The 
effects of altitude and velocity are included in the derivatives with respect to 
those parameters and in the incremental multipliers 

6h = h - h D 

<5V = V - V 0 


where the subscript zero represents the current analysis point (x 0 , x 0 , u Q ) , 
described in the Analysis Point Definition section. 

All stability derivatives are computed as nondimensional terms except the alti- 
tude and velocity parameters. The units of the control derivatives correspond to the 
units used in the nonlinear aerodynamic model. LINEAR is capable of handling both 
Mach and velocity derivatives. Using the equation 

C £m = C £v a h 


(where a^ is the speed of sound at the reference altitude), LINEAR will compute the 
Mach derivatives from the velocity derivatives or vice versa. Internally, however. 
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LINEAR uses velocity derivatives in the computations of the total force and moment 
coefficients. Derivatives with respect to angle of attack and angle of sideslip can 
be obtained in units of reciprocal degrees. These derivatives are simply the cor- 
responding nondimensional derivatives multiplied by 180/ir. 


DATA INPUT 


The interactive linearizer allows the user to input data in various manners : 

1. All input data can be loaded from a single file (either a file assigned to 
logical unit 1 or any file named with 10 characters or less); 

2. Data can be separated onto different files, each containing different sec- 
tions of data; 

3. Data can be typed in from the terminal; or 

4. A combination of (2) and (3). 


Input Files 

The interactive LINEAR input file (defined in table 1) is separated into six 
major sections; vehicle title information; geometry and mass data for the aircraft; 
state, control, and observation variable definitions for the state-space model; the 
trim parameter specification; additional control surfaces that may be specified for 
each case; and various test case specifications. All the input data can be written 
on one file or various files according to the divisions specified at the beginning of 
this paragraph and according to the input format defined in table 1. An example 
input file is listed in appendix F. All the input records to LINEAR are written in 
ASCII form. 


Title Information 

There is one title record that needs to be specified for interactive LINEAR: the 

name of the specific vehicle. This record is read with a 20A4 format. The vehicle 
title appears on each page of the line printer output file. 

Geometry and Mass Data 

The geometry and mass data consist of four records that either can follow the 
vehicle title, exist as a file on its own, or be input interactively from the ter- 
minal. If it is input from the terminal, then it will be stored either on the same 
file as the vehicle title, or on a separate file. The geometry and mass data records 
define the geometry, mass properties, location of the aerodynamic reference point 
with respect to the center of gravity, and angle-of-attack range for the vehicle 
model being analyzed. 
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TABLE 1. - INPUT FORMAT FOR LINEAR 


Input record 





Format 

VEHICLE 

TITLE 





( 20A4 ) 

S b 


c 

Weight 



(4F13.0 ) 

DELX 


z 

DELY 

■^xz 

*xy 

DELZ 

T yz 

LOGCG 


(6F13.0) 
(3F10.0, 12A4) 

“min 

NUMSAT 


“max 

EQUAT 


STAB 



(2F13.0) 

(110, A4, 1 IX, A4) 

STATE 1 


DRVINC 

1 




( 5A4 ,F1 0.0) 

STATE 2 


DRVINC 

2 




( 5A4 , F 1 0 . 0 ) 

STATE 3 

• 


DRVINC 

• 

3 




( 5A4 , F 1 0 . 0 ) 

• 

• 

• 

NUMSUR 


• 

• 





• 

• 

(110) 

CONTROL 

1 

LOCCNT 

1 

CONVR 1 

CNTINC 

1 

(5A4,I10,A4,6X,F10.0) 

CONTROL 

2 

LOCCNT 

2 

CONVR 2 

CNTINC 

2 

( 5 A4 ,110, A4 , 6X , F 1 0 . 0 ) 

CONTROL 

• 

• 

3 

LOCCNT 

• 

« 

3 

CONVR 3 

• 

• 

CNTINC 

• 

• 

3 

( 5A4 ,I10,A4,6X,F10.0) 

• 

• 

• 

NUMYVC 


• 

• 

EQUAT 


• 

• 

• 

• 


• 

• 

(110, A4) 

MEASUREMENT 

1 

PARAM 1 (1-3) 



( 5A4, 3F1 0 . 0 ) 

MEASUREMENT 

2 

PARAM 2 (1-3) 



Ln 

> 

U) 

—A 

O 

• 

O 
*— »* 

MEASUREMENT 

3 

PARAM 3 (1-3) 



(5A4,3F10.0) 


6^ <5„ 6- . 

e min e max a min “max 

NUMXSR 

ADDITIONAL SURFACE 1 
ADDITIONAL SURFACE 2 
ADDITIONAL SURFACE 3 


)y- ■ 5 w- 1 6+-V1 

r min r max cn min tn max 

LOCCNT 1 CONVR1 

LOCCNT 2 CONVR2 

LOCCNT 3 CONVR3 


(8F10.0) 

( IX, 12 ) 

( 5A4 ,110, A4 ) 
(5A4, I 1 0 , A4 ) 
( 5A4 ,110, A4 ) 


ANALYSIS POINT DEFINITION OPTION 
ANALYSIS POINT DEFINITION SUBOPTION 
VARIABLE 1 VALUE 1 

VARIABLE 2 VALUE 2 

VARIABLE 3 VALUE 3 


(20A4) 

(A4) 

( 5A4,F1 5.5) 
( 5A4 ,F1 5.5) 
(5A4,F15.5) 
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TABLE 1 • — Concluded 


Input record 



Format 

NEXT 




(A4) 

ANALYSIS 

POINT 

DEFINITION 

OPTION 

( 20A4 ) 

ANALYSIS 

POINT 

DEFINITION 

SUBOPTION 

( A4 ) 

VARIABLE 

1 

VALUE 1 


(5A4,F15.5) 

VARIABLE 

2 

VALUE 2 


(5A4,F15.5) 

VARIABLE 

3 

VALUE 3 


( 5A4 , FI 5 . 5 ) 


The first and second records describe the vehicle geometry, mass, and mass dis- 
tribution. The first record defines the wing planform area S in units of feet 

squared; the wingspan b in units of feet; the mean aerodynamic chord of the wing c 
in units of feet; and the sea level weight of the vehicle Weight, in units of pounds. 
The second record defines the vehicle moments and products of inertia in units of 
slug-feet squared. The third record defines the offset of the aerodynamic reference 
point with respect to the vehicle center of gravity in the normal right-handed body 
axis reference system with the positive x axis forward. The format of the data can 
be seen in table 1 with DELX, DELY, and DELZ representing the x, y, and z body axis 
displacements of the aerodynamic reference point with respect to the center of gra- 
vity in units of feet (see app. A). The fourth variable of the third record is an 
alphanumeric variable, read using a 12A4 format, to specify if corrections due to a 
center-of-gravity offset are to be computed in LINEAR or in the user-supplied aerody- 
namic subroutine CCALC. The variable LOGCG defaults to a state that causes the aero- 
dynamic reference point offset calculations to be performed by subroutines within 
LINEAR. Any of the following statements in the LOGCG field will cause LINEAR not to 
make these corrections; 

NO CG CORRECTIONS BY LINEAR 

CCALC WILL CALCULATE CG CORRECTIONS 

FORCE AND MOMENT CORRECTIONS CALCULATED IN CCALC 

However, LINEAR will read only the first four characters of the string. 

The final record of this geometry and mass data set defines the angle-of-attack 
range for which the user-supplied nonlinear aerodynamic model (CCALC) is valid. The 
parameters oim^n and a,,,^ specify the minimum and maximum values of angle of attack to 
be used for trimming the aircraft model. These parameters are in units of degrees. 


State, Control, and Observation Variable Definitions 

The state, control, and observation variables to be used in the output formula- 
tion of the linearized system are defined in records that either follow the last of 
the previously described sets of records on the same file, are stored on a separate 
file, or are input through the terminal. If they are input through the terminal, 
they will be stored in a file specified by the user. The number of records in the 
state, control, and observation variable definition set is determined by the number 
of such variables defined by the user. 

The states to be used in the output formulation of the linearized system are 
defined in the first set of records in the state, control, and observation variable 
definitions. The first record of this set, as shown in table 1, defines the number 
of states to be used (NUMSAT), the formulation of the state equation (EQUAT), and 
whether the nondimensional stability derivatives with respect to angle of attack and 
angle of sideslip are to be output in units of reciprocal radians or degrees. The 
variable EQUAT is read using an A4 format and is tested against the following list; 
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NONSTANDARD 

NON-STANDARD 

GENERALIZED 

EXTENDED 

If EQUAT matches the first four characters of any of the listed words, the output 
formulation of the state equation is 

Cx = A ' x + B'u 

If EQUAT is read in as STANDARD or does not match the preceding list, then the de- 
fault standard bilinear formulation of the state equation is assumed, and the out- 
put matrices are consistent with the equation 

X = Ax + Bu 

The variable STAB is also read using an A4 format and is compared with the following 
list: 


DEGREES 

DGR 

If STAB matches the first four characters of either of these words, the nondimen- 
sional stability derivatives with respect to angle of attack and angle of sideslip 
are printed in units of reciprocal degrees on the printer file. Otherwise, these 
derivatives are printed in units of reciprocal radians. 

The remaining records of the state variable definition set are used to specify 
the state variables to be used in the output formulation of the linearized system and 
the increments to be used for the numerical perturbation described in the Linear 
Models section. The state variable names are checked for validity against the state 
variable alphanumeric descriptors listed in appendix D. If a name is not recognized, 
the variable is ignored and a warning message is written to the printer file. The 
increment to be used with any state variable (in calculating the A' and H' matrices) 
and the time derivative of that state variable (in calculating the C and G matrices) 
can be specified using the DRVINC variable. The DRVINC specified for velocity will 
be multiplied by the speed of sound within LINEAR in order to scale up the pertur- 
bation size to a reasonable value while keeping DRVINC on the same order of magnitude 
as for the other states. If DRVINC is not specified by the user, the default value 
of 0.001 is used. 

The next set of records in the state, control, and observation variable defini- 
tions are those defining the variables to be used in the control vector of the output 
model. The first record of this set defines the number of control parameters to be 
used (NUMSUR). The remaining records define the names of these variables (CONTROL), 
their location (LOCCNT) in the common block /CONTRL/ (see the User-Supplied Sub- 
routines section), the units associated with these control variables (CONVR), and 
the increments (CNTINC) to be used with these variables in determining the B' and 
F 1 matrices. 

Because LINEAR has no default control variable names, the control variable names 
input by the user are used for subsequent identification of the control variables. 
Therefore, consistency in the use of control variable names is extremely important. 
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particularly when the user attempts to establish control variable initial conditions 
when using the untrimmed analysis point definition option. 

The CONVR field in the control variable records is used to specify if the control 
variables are given in degrees or radians. CONVR is read using an A4 format and is 
compared to the following list: 


DEGREES 

DGR 

RADIANS 

RAD 

If CONVR agrees with the first four characters of either of the first two names, it 
is assumed that the control variable is specified in units of degrees. If CONVR 
agrees with the first four characters of either of the last two listed names, it is 
assumed that the control variable is specified in units of radians. No units are 
assumed if CONVR does not agree with any of the listed names. When it is assumed 
that the control variable is specified in units of radians, the initial value of the 
control variable is converted to degrees before being written to the printer file. 

The variable CNTINC can be used to specify the increments used for a particular 
control surface when the B' and F' matrices are being calculated. It is assumed that 
the units of CNTINC agree with those used for the surface, and no unit conversion is 
attempted on these increments. If CNTINC is not specified for a particular surface, 
a default value of 0.001 is used. 

The final set of records in the state, control, and observation variable defini- 
tions pertain to the specification of parameters associated with the observation vec- 
tor, observation equation, and observation parameters. The first record of this set 
defines the number of observation variables (NUMYVC) to be used in the output linear 
model and the formulation of the output equation (EQUAT). The remaining records in 
this set specify the variables to be included in the observation vector (MEASUREMENT) 
and any position information (PARAM) that may be required to compute the output model 
for a sensor not located at the vehicle center of gravity. 

The variable used to specify the formulation of the observation equation (EQUAT) 
is compared with the same list of names used to determine the formulation of the 
state equation. If it is determined that the generalized formulation is desired, the 
observation equation 


y=H'x+Gx+F'u 

is used. Otherwise, the standard formulation is assumed, and the form of the obser- 
vation equation used is 


y = Hx + Fu 

The records defining the observation variables to be used in the output formu- 
lation of the linear model contain a variable that includes the parameter name 
(MEASUREMENT) and three fields (PARAM) defining, when appropriate, the location of 
the sensor relative to the vehicle center of gravity. The parameter name is com- 
pared with the list of observation variables given in appendix D. If the parameter 
name is recognized as a valid observation variable name, that observation variable is 
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included in the formulation of the output observation vector. If the parameter name 
is not recognized, an error message is printed and the parameter named is ignored. 
The three variables represented by PARAM(I), PARAM(2), and PARAM( 3 ) provide the x- 
axis, y-axis, and z-axis locations, respectively, of the measurement with respect to 
the vehicle center of gravity if the selected observation is one of the following: 

a x, i 

a y, i 

a z, i 

a n,i 


The unit associated with these variables is feet. If the selected observation var- 
iable is not in the preceding list, the PARAM variables are not used. The sole 
exception to this occurs when Reynolds number is requested as an observation vari- 
able. In that case, PARAM( 1 ) is used to specify the characteristic length. When no 

value is input for PARAM ( 1 ) , the mean aerodynamic chord c is used as the character- 
istic length. 


Trim Parameter Specification 

There is one record in the trim parameter specification set that is associated 
with the subroutine UCNTRL (described in the User-Supplied Subroutines section). 

This record specifies the limits to be used for the trim parameters S e , 6 a , 6 r , and 

6 t h» representing the longitudinal, lateral, directional, and thrust trim parameters, 
respectively. The units associated with these parameters are defined by the imple- 
mentation of UCNTRL. 


Additional Surface Specification 

The first record of the set of additional surface specifications defines the 
number of additional controls to be included in the list of recognized control names 
(NUMXSR). The purpose of defining these additional controls is to allow the user 
to set such variables as landing gear position, wing sweep, or flap position. The 
controls are only defined in the additional surface specification records; actual 
control positions are defined in the analysis point definition records. If such 
controls are defined, the records defining them will have the format specified in 
table 1. The control variable name (ADDITIONAL SURFACE), location (LOCCNT) in the 
common block /CONTRL/, and the units associated with this control variable (CONVR) 
are specified for each additional control. 
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Test Case Specification 


The test case specification records allow the user to define the flight condi- 
tion, or analysis point, at which a linear model is to be derived. Multiple cases 
can be included in the test case specification records. The final record in each set 
directs the program to proceed (NEXT) or to stop (END) execution. 

The first record of a test case specification set determines the analysis point, 
or trim, option to be used for the current case. The ANALYSIS POINT DEFINITION 
OPTION parameter is read in and checked against the list of analysis point defini- 
tion identifiers described in appendix E. The second record of a test case specifi- 
cation set, defining an analysis point definition suboption (ANALYSIS POINT DEFIN- 
ITION SUBOPTION), will be read only if the requested analysis point definition option 
has a suboption associated with it. These suboptions are defined in the Analysis 
Point Definition section. The valid alphanumeric descriptors for these suboptions 
are described in appendix E. 

The remaining records in a test case specification set define test conditions or 
initial conditions for the trimming subroutines. These records consist of a field 
defining a parameter name (VARIABLE) and its initial condition (VALUE). These 
records may be in any order; however, if initial Mach number is to be defined, the 
altitude must be specified before Mach number if the correct initial velocity is to 
be determined. The parameter names are checked against all name lists used within 
LINEAR. Any recognized state, time derivative of state, control, or observation var- 
iable will be accepted. Trim parameters also can be set in these records. 

In general, setting observation variables and time derivatives of the state vari- 
ables has little effect. However, for some of the trim options defined in the Anal- 
ysis Point Definition section, Mach number and load factor are used. The thrust trim 
parameter only affects the specific power trim. For the untrimmed option, the 
initial values of the state and control variables determine the analysis point com- 
pletely. For all other trim options, only certain states are not varied; all con- 
trols connected to the control and engine models are varied. 


Interactive Data Input 

Upon starting program execution, LINEAR will ask the user: 

WHAT IS THE TITLE OF THE CASE BEING RUN? 

The user should answer with some meaningful title followed by a carriage return. 
This title, along with the vehicle title discussed later, will identify this inter- 
active session on all of the output data files. 

The program will then display the following menu: 

PROGRAM LINEAR READS ALL INPUT DATA FROM 1 
UNLESS A DIFFERENT FILE NAME IS SPECIFIED. 


1) 

ALL 

INPUT 

DATA 

IS ON 

1 

2) 

DATA 

IS ON 

ONE 

OTHER 

FILE 
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3) DATA IS ON VARIOUS FILES 

4) INPUT DATA FROM A FILE(S) AND 
FROM THE TERMINAL 

5) INPUT ALL DATA FROM THE TERMINAL 

The user must answer with the appropriate number followed by a carriage return. 

If the user answers with a "1", the program will read all the input data from the 
file assigned logical unit one. If the user selects choice "2", the program will ask 
for the file name that contains all the input data required by LINEAR. If the user 
selects "3" or "4", the program will prompt the user for each file name pertaining to 
each section of input data. If option "5" is chosen, the program will prompt the 
user for the input data values. 

The data is separated into the six sections described previously: 

1 . Vehicle title information, 

2. Vehicle geometry and mass properties, 

3. State, control, and observation variable definitions, 

4. Limits for the trim parameters, 

5. Additional control surface definitions, and 

6. Test case specifications. 

The program will prompt the user for the name of the file that contains each section 
of data. An example prompt follows: 

WHAT IS THE NAME OF THE FILE CONTAINING THE TITLE INFORMATION? 

(TYPE 0 TO INPUT THE TITLES FROM THE TERMINAL AND STORE ON A FILE) 

If the user already has created a file containing the appropriate data needed for the 
section, a 10-character (maximum) file name may be entered. If the user wishes to 
type the data from the terminal and have the program format a file, a zero may be 
entered and the program will prompt 

WHAT WOULD YOU LIKE THE FILE TO BE CALLED? 

At this time the user must input a file name (10 characters maximum) on which all 
data typed from the terminal will be stored for that section of input data. LINEAR 
will then prompt the user for the specific data it needs for this section. 

If the user wishes to input all the data from the terminal, selection "5" should 
be chosen and the program will prompt for a data file name on which to store the 
input data. The program will prompt the user for all the input data required. 
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When inputting the title data, the program reads character strings of 80 char- 
acters. When prompting for vehicle mass and geometry data, the program reads a 
floating point field. An example data prompt is as follows: 

INPUT THE VEHICLE GEOMETRY AND MASS PROPERTIES 
WING AREA (FT**2) = 

The user would input the vehicle wing area in square feet followed by a carriage 
return. The remainder of the vehicle geometry data is input similarly. 

When inputting the state vector from the terminal, the user must define the size 

and formulation of the state equation as well as the units in which the angle of 

attack and sideslip stability derivatives will be output. The program will prompt 

HOW MANY STATE VARIABLES ARE IN YOUR MODEL? 

(MAXIMUM 12) 

The program will read an integer variable and then prompt 

WHICH FORM OF THE STATE EQUATION DO YOU 
WISH TO USE? 

1 ) DX/DT = A*X + B*U + D*V 

2) C * DX/DT = A*X + B*U + D*V 

If the standard formulation is being used, the user should input "1" and for the 
generalized equation, "2". The LINEAR program uses the generalized formulation for 
internal calculations and then performs any transformations necessary for the output 
formulation requested. The program will then prompt for the units of the angle of 
attack and sideslip stability derivatives. The program uses radians for all inter- 
nal calculations and transforms the stability derivatives to degrees for output 
if desired. 

The variables used in the state vector are input next and are read in strings 
20 characters in length. Each variable name entered is checked against legal names 
for state variables. If a variable is incorrect, the program will respond 

VARIABLE so-and-so 

IS AN INVALID STATE PARAMETER 

PLEASE CHOOSE ANOTHER OR TYPE HELP FOR MORE INFO. 

If the user types HELP, the program will list all the valid names for the state 
variables. The names are given in appendix D. 

After the entire state vector has been entered, the program will prompt 

THE STEP SIZE USED FOR DERIVATIVE EXTRACTION 
IS INITIALIZED FOR ALL VARIABLES TO 0.001. 

DO YOU WISH TO INPUT A DIFFERENT VALUE 
FOR ANY OF THE STATE VARIABLES? (Y/N) 
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To change the perturbation step size for any of the state variables, the user should 
enter "Y" and, after being prompted, the variable name to be updated and the desired 
step size. 

For the interactive input of the control variables, similar prompts are given 
with one major difference: after each variable name is input, the program prompts 

for the location of the variable in the common block /CONTRL/ (the named common block 
/CONTRL/ reserves 30 locations for control surface variable values used as interface 
to the user-supplied aerodynamic model). The program also prompts for the units to 
be used in calculations involving the control surface: 

WHAT UNITS ARE USED IN CALCULATIONS WITH 
THIS CONTROL SURFACE? 


1) 

NONE 

2) 

DEGREES 

3) 

RADIANS 


The units should be consistent with the calculations in UCNTRL. The program performs 
transformations based on the surface type for the printer file only. No other con- 
trol surface transformations are performed. Again, the derivative step size can be 
changed for individual variables specified by the user. 

For the interactive input of the observation vector, similar prompts are given 
with a few variations. The forms available for the observation equation are 

1 ) Y = H*X + F*U + E*V 

2) Y = H*X + G*DX/DT + F*U + E*V 

For variables defined in the observation vector as being offset from the vehicle cen- 
ter of gravity, the program will prompt for the offset in the x, y, and z stability 
axes. For any variables that are calculated based on some reference value, such as 
Reynolds number, the program prompts for the correct value when the name is entered. 
Again, the step size for derivative extraction can be modified for any observation 
variable. Valid observation variable names are listed in appendix C. 

If selected to be input from the terminal, the program will prompt the user for 
the limits to be used for the trim parameters 6 e , 6 a , 6 r , 6^h' which represent the 
longitudinal, lateral, directional, and thrust trim parameters, respectively. The 
units associated with these parameters are defined by the implementation of UCNTRL. 
(See app. J of this report for details.) An example prompt is as follows: 

WHAT IS THE MAXIMUM THRUST PARAMETER? 

Again, the program reads a floating point field. 

After the trim axis parameter limits are defined, the program prompts 

HOW MANY ADDITIONAL SURFACES DO YOU WISH 
TO DEFINE? 
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At this time, the user can define additional controls to be included in the list of 
recognized control names. The purpose of defining these additional surfaces is to 
allow the user to set such variables as landing gear position, wing sweep, or flap 
position. The control surface name, location in common /CONTRL/, and units are 
defined here for each surface; actual control surface positions are defined during 
analysis point definition. 

The program will then prompt the user to indicate what aerodynamic model will 
be used for this analysis. Aerodynamic models are described in the next section of 
this report. 

If a nonlinear aerodynamic model is used, the program will prompt the user to 
determine if there are additional control surfaces for which derivatives are desired. 
These controls will not become part of the control vector. This option of inter- 
active LINEAR does not allow these control variables to be read from or saved to a 
file, and is not available on the batch version of LINEAR. However, the stability 
and control derivative output file will contain these surfaces and their derivatives. 

The point at which the nonlinear system equations are linearized is referred to 
as the analysis or trim point. This can represent a true steady-state condition on 
the specific trajectory (a point at which the rotational and translational accelera- 
tions are zero) or a totally arbitrary state on a trajectory. LINEAR allows the user 
to select from a variety of analysis points. Within the program, these analysis 
points are referred to as trim conditions, and several options are available to 
the user. 

If the user is defining the trim conditions interactively, instead of reading the 
data from a formatted file, the program will prompt: 

WHAT TYPE OF TRIM DO YOU WISH TO RUN? 

1 ) NO TRIM 

2) STRAIGHT AND LEVEL FLIGHT 

3 ) LEVEL TURN 

4 ) PUSH-OVER/PULL-UP 

5) THRUST STABILIZED TURN 

6 ) BETA 

7) SPECIFIC POWER 

The user must specify the number corresponding to the type of trim desired. The 
arbitrary state and control option is designated as NO TRIM, and in selecting this 
option the user must specify all nonzero state and control variables. 

For all of the analysis point definition options, any state or control parameter 
may be input after the required data is defined when the program prompts as follows: 

WHAT OTHER STATE OR CONTROL VARIABLES WOULD 
YOU LIKE TO INITIALIZE? 

TYPE HELP TO LIST VALID VARIABLES NAMES 

(ENTER "N" IF NO OTHER CONDITION IS TO BE SPECIFIED) 
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The user should enter the name of the variable to be initialized or type HELP to 
list the variables that may be initialized. 

If a level turn trim is selected, the user can specify the direction of the turn 
by entering "TURN RIGHT" or "TURN LEFT" when prompted to initialize any other var- 
iables. This is especially useful for asymmetric aircraft. Default is "TURN RIGHT". 

If data are being read from a file(s), the data can be reviewed by replying "Y" 
when prompted 

DO YOU WISH TO REVIEW THE DATA AS IT 
IS READ FROM THE INPUT FILE(S)? (Y/N) 

If any data are entered incorrectly, the user should finish entering all data as 
prompted by the program until asked if there are any changes to be made: 

DO YOU WISH TO CHANGE ANY PARAMETERS 
IN YOUR MODEL OR UPDATE THE INPUT 
CASE SELECTION FILE? (Y/N) 

At this time, the user selects "Y" and the program prompts 

WHAT PARAMETERS DO YOU WISH TO CHANGE? 

1 ) VEHICLE GEOMETRY 

2 ) STATE VECTOR 

3) OBSERVATION VECTOR 

4) CONTROL VECTOR 

5) CASE INPUT FILE 

6) STABILITY AND CONTROL DERIVATIVES 

7) SURFACES TO BE SET 

8) NO FURTHER CHANGES 

Any section of data can be updated by reading a new data file or by typing in the 
data from the terminal. The program prompts the user to determine the source of 
input as explained previously. If the user wishes to modify the vehicle geometry and 
mass data or the stability and control derivatives interactively, the program saves 
all the values from the previous case and allows the user to modify any individual 
parameter. For example, if the user wishes to modify the vehicle geometry proper- 
ties, a "1" would be entered and the program would continue 

WHAT IS THE NAME OF THE FILE CONTAINING THE VEHICLE 
GEOMETRY AND MASS PROPERTIES? 

(TYPE 0 TO INPUT THE DATA FROM THE 
TERMINAL AND STORE ON A FILE) 

If the user replies with a zero, the program will prompt for a new output file name 
and then respond as follows: 
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GEOMETRY AND MASS DATA FOR: 
USER’S GUIDE 

FOR THE PROJECT: 

DEMONSTRATION CASE 


1 ) 

WING AREA 

608.000 

2) 

WING SPAN 

42.800 

3) 

MEAN CHORD 

15.950 

4) 

VEHICLE WEIGHT 

40700.000 

5) 

IX 

28700.000 

6) 

IY 

165100.000 

7) 

IZ 

187900.000 

8) 

IXZ 

-520.000 

9) 

IXY 

0.000 

10) 

IYZ 

0.000 


VECTOR DEFINING REFERENCE POINT OF AERODYNAMIC MODEL 
WITH RESPECT TO THE VEHICLE CENTER OF GRAVITY: 


1 1 ) DELTA X 0.000 

12) DELTA Y 0.000 

13) DELTA Z 0.000 

14) MINIMUM ANGLE OF ATTACK -10.000 

15) MAXIMUM ANGLE OF ATTACK 40.000 


16) FORCE AND MOMENT COEFFICIENTS CORRECTIONS DUE 
TO THE OFFSET OF THE REFERENCE POINT OF THE 
AERODYNAMIC MODEL FROM THE VEHICLE CG ARE 
INTERNALLY CALCULATED IN LINEAR. TYPE IN 
16=0 IF CALCULATIONS ARE DONE IN CCALC. 


TO MAKE CHANGES, TYPE IN THE LINE NUMBER ON THE 
DESIRED VALUES. 

"N" INDICATES NO CHANGES ARE DESIRED 
"R" REFRESHES THE SCREEN WITH NEW VALUES. 

FOR EXAMPLE: TO CHANGE VEHICLE WEIGHT 

4=40000.0 

If selection "16=0" is chosen, the program will set a logical variable that deter- 
mines where center-of- gravity corrections are made to the force and moment coef- 
ficients. The vector defining the reference point of aerodynamic model with respect 
to the vehicle center of gravity must be defined in feet. 

If the state, control, or observation vectors need to be modified, the user must 
input the entire vector again. The choices available for the test case selection are 
to input the data from the terminal for each case or read the entire case specifica- 
tion from a file. 
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All surfaces in the control model are initially zero unless specified in the case 
input file. Surfaces not in the control vector can be set by choosing option "7" 
when prompted for parameter changes. The program will prompt for the file name con- 
taining the surfaces that can be set for each case. If the user wishes to input the 
surface names interactively, the program will prompt for the number of surfaces to be 
added, their names, location in common /CONTRL/ and units used in subroutine UCNTRL. 
If the location in common /CONTRL/ has already been designated for another surface, 
the program will save the original name and ignore any other name input for that 
location, and so inform the user. 

The interactive LINEAR program provides the user with two different ways to 
calculate the total force and moment coefficients. One method is by using a full 
aerodynamic model and the subroutines ADATIN and CCALC. The other method uses a set 
of linearized stability and control derivatives, and the internal subroutine CLNCAL. 
When using the second method, the data can be typed in during execution or can be 
input from a file. 

An example of a linearized stability and control data set is shown in appendix G. 
The first line contains the file description. The case number (1-999) at the end of 
the first line allows the user to correlate the linearized data from a case to other 
data from that case. The next line is the title of the run. The following two lines 
of text describe the project for which the case was run, and the next two lines list 
the reference altitude and Mach number for the case. The next line determines 
whether the alpha and beta derivatives are in units of degrees or radians. The 
following two lines indicate whether Mach or velocity derivatives will be read in. 
This is intended to warn the user against modifying the wrong derivative while 
editing outside of LINEAR. If the wrong derivative is changed, LINEAR will ignore 
the change. Additionally, if the reference altitude is changed, LINEAR will rely 
on the specified derivative and compute the other based on the speed of sound which 
is a function of altitude. While editing this file in LINEAR, this caution is not 
necessary since LINEAR will perform the computation right away. Also, when these 
derivatives are changed in LINEAR, the last velocity or Mach derivative changed will 
determine which one LINEAR will regard as the specified derivative when the file 
is output. 

The next line determines the number of control surfaces used for control deriva- 
tives. The units of each of the control surfaces, where applicable, are listed next 
to the surface name as well as the location of the surface in the common /CONTRL/. 

The units defined for each surface correspond to the units defined in user-supplied 
subroutine CCALC. 

The rest of the file contains the stability and control derivatives. This part 
is broken into six sections, one for each of the force and moment coefficients. In 
each section, there is one zero coefficient corresponding to zero angle of attack, 
zero sideslip, zero control surface deflection, and constant Mach and altitude. Each 
section also contains 10 stability derivatives. In addition, each section has as 
many control derivatives as specified earlier. The six sections must always be in 
the same order: rolling moment, pitching moment, yawing moment, drag force, lift 

force, and sideforce. 

The easiest way to format a derivative file is by running the interactive LINEAR 
program and inputting the derivatives when prompted by the program. One of the out- 
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put files from LINEAR is the stability and control derivative file that can be used 
later to analyze the same case. (See the Aerodynamic Model section for more infor- 
mation. ) When using a linearized set of stability and control derivatives from a 
previous case, the user must be careful to define a control vector that is consistent 
with the derivative set. If a surface is defined that does not have corresponding 
control derivatives on the file, that surface's effects on the forces and moments 
will not be included in the linear system model calculated. The user should also be 
aware that LINEAR expects to have control over all axes to determine a trimmed state. 
If there are no nonzero control derivatives in a given axis, it is highly likely that 
LINEAR will fail to trim; and it may even generate a fatal error in the nonlinear 
equation solver (generally a floating-point overflow). 


AERODYNAMIC MODEL 


As already stated, there are two methods for calculating the total force and 
moment coefficients. One method relies on a full nonlinear aerodynamic model to pro- 
vide the stability and control derivatives. These derivatives are stored in common 
blocks in the user-supplied routine ADATIN and used in CCALC to calculate the total 
force and moment coefficients using the state and control surface values calculated 
in LINEAR. These subroutines and the aerodynamic model are the same ones used in the 
manned simulations at Ames-Dryden. 

Using a full aerodynamic model with CCALC and ADATIN enables the user to run as 
many cases as desired without changing derivative data sets. From each case, the 
program generates a structured linearized derivative data file which can be used as 
an input file later. This file contains all the stability and control derivatives 
for the control surfaces specified in the control vector, and any other controls 
specified by the user. The interface required for using a nonlinear aerodynamic 
model is described later in the Aerodynamic Model Subroutines section. 

The other method for calculating total force and moment coefficients requires the 
user to supply linearized stability and control derivatives. These are used in the 
internal subroutine CLNCAL to calculate the total force and moment coefficients. 

When using this option, care must be taken to ensure that the input data is valid 
for the trim point desired. Center-of-gravity shifts will not affect the data that 
is output to the derivative file when this method is used. 

The method used to calculate the force and moment coefficients is specified 
during the input phase for a particular case. After all the vehicle data, state, 
control, observation vectors, and case data have been defined, the program will 
prompt 


DO YOU WISH TO 

1 ) INPUT A SET OF LINEAR STABILITY 
AND CONTROL DERIVATIVES 

2) OR USE THE NONLINEAR MODEL 

YOU HAVE INTERFACED WITH THIS PROGRAM 


If the internal routines are to be used to calculate the derivatives, the program 
will prompt 
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DO YOU WISH TO 


1 ) GENERATE A NEW DATA SET 

2 ) OR USE A DERIVATIVE DATA SET YOU HAVE ON FILE 

3) OR MODIFY THE OLD ONE 

J 

If the derivatives are stored on a file (format for the file is given in app. G), the 
program will prompt for a file name. At this point, the user can modify the data 
interactively by choosing one of the following options: 

TO CHANGE ANY DERIVATIVES, INPUT THE 

APPROPRIATE LINE # TO SEE THE DATA. 

TYPE IN "N" IF NO CHANGES DESIRED 

1 ) TYPE OF ANGLE MEASUREMENT 

2 ) ROLL MOMENT 

3) PITCH MOMENT 

4) YAW MOMENT 

5) DRAG FORCE 

6 ) LIFT FORCE 

7) SIDE FORCE 

8 ) ADD SURFACES 

If the user selects option one, the program will list the units that are used for 
output of the stability derivatives. It will also list the units for each control 
surface derivative as used for calculations in the subroutine UCNTRL . There are 
three choices for units for the control surfaces: none, radians, and degrees. If 

the control surface calculations in UCNTRL are nondimensional, the user should spec- 
ify none, and no conversion is done. For surface calculations in radians, a conver- 
sion to degrees for output only is performed. If degrees is specified, no conver- 
sion is done. When the user inputs the stability and control derivatives, the pro- 
gram does not perform any conversions, therefore, care must be taken to ensure that 
the units are consistent with the aerodynamic model and the calculations in UCNTRL. 

If option "8" is chosen, the control surfaces defined are listed, and the user 
can define other surfaces in common /CONTRL/ for which derivatives exist. 

If the derivatives are be to calculated from a nonlinear aerodynamic model, the 
program will prompt for the names of any other control surfaces not defined in the 
control vector for which linearized derivatives will be calculated. If options "2" 
through "7" are chosen, the program will list the current values for the stability 
and control derivatives for each respective coefficient. The user can change these 
values interactively, as in the following example: 

ROLL MOMENT DERIVATIVES 


WITH RESPECT TO: 


0) 

CLO 


-0.000033820 

1 ) 

ROLL RATE 

(RAD/SEC) 

-0.200000003 

2) 

PITCH RATE 

(RAD/SEC) 

0.000000000 

3) 

YAW RATE 

(RAD /SEC) 

0.163589999 
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4) 

VELOCITY 

(FT/ SEC) 

-0.000000115 

5) 

MACH NUMBER 


-0.000119522 

6) 

ALPHA 

(DEG/SEC) 

0.000362335 

7) 

BETA 

(DEG/SEC) 

-0.140983000 


• • 


TO MAKE CHANGES, TYPE IN THE LINE NUMBER 
AND THE VALUE OF THE DERIVATIVE. 

"N" INDICATES NO CHANGES ARE DESIRED 

"R" REFRESHES THE SCREEN WITH THE NEW VALUES 

FOR EXAMPLE: ROLL MOMENT WRT ALPHA 

6=0.003 


OUTPUT FILES 


There are four output files from LINEAR: (1) a general-purpose analysis file, 

(2) a printer file containing the calculated case conditions and the state-space 
matrices for each case, (3) a printer file containing the calculated case condition! 
only, and (4) a linearized set of stability and control derivatives. LINEAR also 
outputs any data input by the user during an interactive session that can be used 
later either by the batch or interactive programs . The program prompts the user fo. 
the file names that each of the sections of data will be written to, as described i 
the Interactive Data Input section. 

The general-purpose analysis file contains the following: the title of the cas 

being analyzed; the state, control, and observation variables used to define the 
state-space model; and the state and observation matrices calculated in LINEAR. Th 
C and G matrices are printed only if the user has selected an appropriate formula- 
tion of the state and observation equations. The output for this file is assigned 
to FORTRAN device number 15. An example general-purpose analysis file is presented 
in appendix H corresponding to the format of table 2. 

The vehicle and case titles are written on the first two records of the analysi 
file in 80-character strings and are specified in LINEAR as the title of the vehicl 
and the title of the cases. The next record contains the number of the case as 
defined in LINEAR (maximum of 999 cases). 

The number of states, controls, and outputs used to define each case is written 
on the subsequent record. The formulation of the state and observation equations 
is listed next, followed by the names and values of the states, controls, dynamic 
interaction variables, and outputs. These values are followed by the matrices that 
describe each case. 

The title records only appear at the beginning of the file while all other 
records are duplicated for each subsequent case calculated in LINEAR. The matrices 
are written row-wise, five columns at a time as shown in table 3. This table shows 
system containing seven states, three controls, and eleven outputs using the genera 
state equation and standard observation equation. 
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TABLE 2. - ANALYSIS FILE FORMAT 


Variable 

Format 


Title of the case 

(4A20) 


Title of the aircraft 

(4A20) 


Case number 

(// ,64X,I3 ) 


Number of states, controls, and outputs 

(17X,I2,22X,I2,22X, 

13) 

State equation formulation 

( 36X, 2A4 ) 


Observation equation formulation 

( 36X, 2A4 ) 
(////) 


State variable names, values and units 

( IX, 5A4 ,3X,E12.6,2X 
(////) 

, A20 ) 

Control variable names, values and units 

( 1X,5A4,3X,E12.6,2X 
(////) 

, A20 ) 

Dynamic interaction variable names and units 

(1X,5A4,17X,A20) 

(////) 


Output variable names, values and units 

( IX, 5A4 ,3X,E12.6,2X 

, A20 ) 

Matrix name 

( // , A8 , / ) 


A-Matrix 

(5 (E13.6) ) 


Matrix name 

(// ,A8,/) 


B-Matrix 

( 5 (E 1 3 . 6 ) ) 


Matrix name 

(// ,A8,/) 


D-Matrix 

( 5 (El 3 . 6 ) ) 


Matrix name 

( // , A8 , / ) 


C-Matrix (if general form chosen) 

( 5 (El 3 . 6 ) ) 


Matrix name 

(// , A8 ,/ ) 


H-Matrix 

(5 (E13 .6 ) ) 


Matrix name 

( // , A8 , / ) 


F-Matrix 

(5 (E13 .6 ) ) 


Matrix name 

( // , A8 , / ) 


E-Matrix 

(5 (E13.6 ) ) 


Matrix name 

(// ,A8,/) 


G-Matrix (if general form chosen) 

(5(E13 .6 ) ) 




TABLE 3. - ANALYSIS FILE OUTPUT 
MATRIX STRUCTURE 


Size 

of i 

matrix 

Output 

formulation 

A = 

[ 7 

X 

7] 

A 

= 

[ 7 

X 

5] 







[ 7 

X 

2] 

B = 

[ 7 

X 

3] 

B 

= 

[ 7 

X 

3] 

D = 

t 7 

X 

6] 

D 

= 

[ 7 

X 

5] 







[ 7 

X 

1] 

C = 

[ 7 

X 

7] 

C 

= 

[ 7 

X 

5] 







[ 7 

X 

2] 

H = 

[11 

X 

7] 

H 

= 

[11 

X 

5] 







[11 

X 

2] 

F = 

[11 

X 

3] 

F 

= 

[11 

X 

3] 

E = 

[11 

X 

6] 

E 

= 

[11 

X 

5] 







[11 

X 

1] 
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The second output file, which is assigned to FORTRAN device unit number 3, con- 
tains the values calculated in LINEAR describing each case. The first section of 
I this file contains a listing of the input data defining the aircraft's geometry and 
mass properties, variable names defining the state-space model, and various control 
surface limits characteristic of the given aircraft. Appendix I presents an example 
printer output file. " 

The second section of this file contains the trim conditions of the aircraft at 
the point of interest. These conditions include the type of trim being attempted, 

whether trim was achieved, and parameters defining the trim condition. 

i 

The values for the variables of the state-space model at the trim condition are 

• • • • ♦ • 

• also printed. If trim was not achieved, p, q, r, V, 3, and a (calculated from the 
' equations of motion) and the force and moment coefficients are printed. Changes in 
the geometry and mass properties are also printed. 

The third section of this output file contains the nondimensional stability and 
control derivatives for the trim condition calculated. The static margin of the 
aircraft at the given flight condition is also printed. 

1 The final section of this output file contains the state and observation matrices 

for the given flight condition. The formulation of the state equations and only the 
terms of the matrices chosen by the user to define the model are printed. 

The third output file, which is assigned to FORTRAN device unit 2, contains the 
trim conditions of the aircraft at the point of interest. These conditions include 
the type of trim being attempted, whether trim was achieved, parameters defining the 
trim condition, and the static margin of the aircraft at the given flight condition. 
Appendix I also presents an example of this file. 

The fourth output file is the set of linearized stability and control derivatives 
for each case as defined in the Aerodynamic Model section. 


USER-SUPPLIED SUBROUTINES 


There are five subroutines that must be supplied by the user to interface LINEAR 
with a specific aircraft's subsystem models: ADATIN, CCALC, IFENGN, UNCTRL, and 

MASGEO. The first two subroutines comprise the aerodynamic model. The subroutine 
IFENGN is used to provide an interface between LINEAR and any engine modeling rou- 
tines the user may wish to incorporate. UCNTRL converts the trim parameters used 
by LINEAR into control surface deflections for the aerodynamic modeling routines. 
MASGEO enables the mass and geometry properties to change as the aircraft configura- 
tion does (for example, with wing sweep). The use of these subroutines is illustra- 
ted in figure 3 which shows the program flow and the interaction of LINEAR with the 
user- supplied subroutines. These subroutines are described in detail on the fol- 
lowing pages. Examples of these subroutines can be found in appendix J. 


53 



Figure 3. Program flow diagram 
showing communication with user- 
supplied subroutines. 


Aerodynamic Model Subroutines 

It is assumed that the aerodynamic models consist of two main subroutines: 

ADATIN and CCALC. ADATIN is used to input the basic aerodynamic data from remote 
storage. ADATIN can also be used to define aerodynamic data. CCALC is the sub- 
routine that uses this aerodynamic data, the state variables, and the surface posi- 
tions to determine the aerodynamic coefficients. Either routine may call other 
subroutines to perform related or required functions; however, from the point of 
view of the interface to LINEAR, only these two subroutines are required for an aerc 
dynamic model. 

Externally, ADATIN has no interface to the program LINEAR. The subroutine is 
called only once when the aerodynamic data is input or defined. The calling progran 
has to provide ADATIN with the input devices it requires, but no other accommodation 
is necessary. The aerodynamic data is communicated from ADATIN to CCALC through 
named common blocks that occur only in these two routines. 

The interface between CCALC and the calling program is somewhat more complicated 
However, the interface is standard and this feature provides a framework about which 
a general purpose tool can be built. This interface consists of several named com- 
mon blocks which are used to pass state variables, air data parameters, surface posi 
tions, and force and moment coefficients between CCALC and the calling program. 

CCALC is executed whenever new aerodynamic coefficients are required (for example, 
once per frame for a real-time simulation). 

The main transfer of data into the subroutine CCALC is through five named common 
blocks. These common blocks contain the state variables, air data parameters, and 
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Externally, ADATIN has no interface to the program LINEAR. The subroutine is 
called only once when the aerodynamic data is input or defined. The calling program 
has to provide ADATIN with the input devices it requires, but no other accommodation 
is necessary. The aerodynamic data is communicated from ADATIN to CCALC through 
named common blocks that occur only in these two routines. 

The interface between CCALC and the calling program is somewhat more complicated. 
However, the interface is standard and this feature provides a framework about which 
a general purpose tool can be built. This interface consists of several named com- 
mon blocks which are used to pass state variables, air data parameters, surface posi- 
tions, and force and moment coefficients between CCALC and the calling program. 

CCALC is executed whenever new aerodynamic coefficients are required (for example, 
once per frame for a real-time simulation) . 

The main transfer of data into the subroutine CCALC is through five named common 
blocks. These common blocks contain the state variables, air data parameters, and 
surface positions. The transfer of data from CCALC is through a named common block 
containing the aerodynamic force and moment coefficients. The details of these com- 
mon blocks follow. 


The common block /DRVOUT/ contains the state variables and their derivatives with 
respect to time. The structure of this common block is given as follows: 

COMMON /DRVOUT/ T , 


p 


Q 

R 

V 


ALP , 

BTA 

THA 


PSI , 

PHI 

H 


X 

Y 

TDOT 




PDOT 


QDOT , 

ROOT 

VDOT 


ALPDOT , 

BTADOT 

THADOT 


PS I DOT, 

PHI DOT 

HDOT 


XDOT , 

YDOT 


The body axis rates p, q, and r appear as P, Q, and R, respectively. Total velocity 
is represented by the variable V and the altitude by H. Angle of attack (ALP), angle 
of sideslip ( BTA ) , and their derivatives with respect to time (ALPDOT and BTADOT, 
respectively) are also contained within this common block. 


The common block /SIMOUT/ contains the main air-data parameters required for the 
aerodynamic model. The variables in this common block are as follows: 


COMMON /SIMOUT/ AMCH, 

QBAR, 

GMA 

UB , 

VB , 

WB 

VEAS, 

VC AS 



Mach number and dynamic pressure are the first two entries in the common block sym- 
bolized by AMCH and QBAR, respectively. The body axis velocities u, v, and w are 
included as UB, VB, and WB, respectively. 


55 



The /CONTRL/ common block contains the surface position information. The exact 
definition of each of the elements used for a particular aerodynamic model is deter 
mined by the implementer of that model. The structure of the common block /CONTRL/ 
is as follows: 


COMMON /CONTRL/ DC (30) 

The common block /DATAIN/ contains geometry and mass information. The first 
three variables in the common block, S, B, and CBAR, represent wing area, wing span 
and mean aerodynamic chord, respectively. The vehicle mass is represented by AMSS. 
The structure of common /DATAIN/ is as follows: 


COMMON /DATAIN/ S , 

B 

CBAR, 

AMSS 

AIX , 

AIY , 

AIZ , 

AIXZ 

AIXY , 

AIYZ , 

AIXE 



The information on the displacement of the reference point of the aerodynamic 
data with respect to the aircraft center of gravity is contained in the /CGSHFT/ 
common block: 


COMMON /CGSHFT/ DELX, DELY, DELZ 

The displacements are defined along the vehicle body axis with DELX, DELY, and DELZ 
representing the displacements for the x, y, and z axes, respectively. 

The output common block /CLCOUT/ contains the variables representing the aero- 
dynamic moment and force coefficients: 

COMMON /CLCOUT/ CL, CM, CN, CD, CLFT, CY 

The variables CL, CM, and CN are the symbols for the rolling moment (C^), pitching 

moment (C m ), and yawing moment (C n ) coefficients, respectively; these terms are bod 
axis coefficients. The stability axis forces are represented by CD, coefficient of 
drag (Cq); CLFT, coefficient of lift (C L ); and CY, sideforce coefficient (C Y )* 


Control Model Subroutines 

The LINEAR program attempts to trim the given condition using control inputs 
similar to that of a pilot: longitudinal stick, lateral stick, rudder, and throttl 

The UCNTRL subroutine uses these trim output control values to calculate actual sur 
face deflections and power lever angles for the given aircraft (fig. 2). The loca- 
tion of each surface and power lever angle in the /CONTRL/ common block is specific 
by the user in the input file (maximum of 30 surfaces). The limits for the control 
parameters in pitch, roll, yaw, and thrust are user specified and must agree in uni 
with the calculations in CCALC. 

The common block /CTPARM/ contains the four trim parameters that must be relate 
to surface deflections in the subroutine UCNTRL: 

COMMON /CTPARM/ DES , DAS, DRS , THRSTX 
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I The output from UCNTRL is via the common block /CONTRL/ described previously in the 
Aerodynamic Model Subroutines section. The variables DES, DAS, DRS, and THRSTX are 
used to trim the longitudinal, lateral, directional, and thrust axes, respectively. 

For an aircraft using feedback (such as a statically unstable vehicle), state 
variables and their derivatives are available in the common block /DRVOUT/. However, 
this control model must contain no states of its own. If parameters other than state 
variables and their time derivatives are required for feedback, the user may access 
them using the common block /OBSERV/ described later in the Mass and Geometry Model 
Subroutines section of this report. 


Engine Model Subroutines 

, The subroutine IFENGN computes individual engine parameters used by LINEAR to 
calculate force, torque, and gyroscopic effects due to engine offset from the cen- 
terline. The parameters for each engine (maximum of four engines) are passed through 
common /ENGSTF/ as follows: 

COMMON /ENGSTF/ THRUST (4), TLOCAT(4,3), XYANGL(4 ) , XZANGL(4), 

TVANXY ( 4 ) , TVANXZ ( 4 ) , DXTHRS ( 4 ) , 

EIX (4), AMSENG (4), ENGOMG ( 4 ) 

The variables in this common block correspond to thrust (THRUST); x, y, and z body 
axis coordinates of the point at which thrust acts (TLOCAT); the orientation of the 
thrust vector in the x-y body axis plane (XYANGL) , in degrees; the orientation of the 

thrust vector in the x-z body axis plane (XZANGL), in degrees; the orientation of the 

thrust vector in the x-y engine axis plane (TVANXY), in degrees; and the orientation 

of the thrust vector in the x-z engine axis plane (TVANXZ), in degrees. The distance 

between the center of gravity of the engine and the thrust point (DXTHRS), measured 
positive in the negative x engine axis direction; the rotational inertia of the 
engine (EIX); mass (AMSENG); and the rotational velocity of the engine (ENGOMG) 

'are also in this common block. The variables are all dimensioned to accommodate up 
to four engines. Appendix B provides graphical and physical descriptions of these 
parameters . 

While the common block structure within LINEAR is designed for engines that do 
not interact with the vehicle aerodynamics, propeller-driven aircraft can be easily 
modeled by communicating the appropriate parameters from the engine model IFENGN to 
the aerodynamic model CCALC. 


Mass and Geometry Model Subroutines 

The subroutine MASGEO allows the user to vary the center-of-gravity position and 
vehicle inertias automatically. Nominally, this subroutine must exist as one of the 
user subroutines, but it may be nothing more than a RETURN and END statement. MASGEO 
is primarily for variable- geometry aircraft (such as an oblique-wing or variable- 
sweep configuration) or for modeling aircraft that undergo significant mass or iner- 
tia changes over their operating range. The center-of-gravity position and inertias 
may be functions of flight condition or any surface defined in the /CONTRL/ common 


57 


block. Changes in center-of-gravity position are passed in the /CGSHFT/ common 
block, and inertia changes are passed in the /DATAIN/ common block. 

Care must be taken when using the subroutine MASGEO in combination with select! 
an observation vector that contains measurements away from the center of gravity, 
measurements are desired at a fixed location on the vehicle, such as a sensor plat- 
form or nose boom, the moment arm to the new center-of-gravity location must be 
recomputed as a result of the center-of-gravity shift for accurate results. This c 
be accomplished by implementing the moment arm calculations in one of the user 
subroutines and passing the new moment arm values through the /OBSERV/ common block 

COMMON /OBSERV/ OBVEC (120), PARAM (120,6) 

The common block /OBSERV/ allows the user to access all the observation vari- 
ables during trim as well as to pass parameters associated with the observations 
back to LINEAR. The common block /OBSERV/ contains two single- precision vectors: 
OBVEC(120), and PARAM( 1 20 , 6 ) . A list of the available observations and parameters 
given in appendix C. Access to the observation variables allows the users to imple 
ment trim strategies that are functions of observations, such as gain schedules anc 
surface management schemes. The parameters associated with the observations are us 
primarily to define the moment arm from the center of gravity to the point at which 

the observation is to be made. If the user subroutine MASGEO is used to change the 

center-of-gravity location and observations are desired at fixed locations on the 
vehicle, then the moment arm from the new center-of-gravity location to the fixed 
sensor location ((x, y, z), in feet) must be computed in one of the user subroutine 
and passed back in the first three elements of the PARAM vector associated with the 
desired observation ( PARAM ( n = 1 to 3), where n is the number of the desired observ 

tion). Additional information on observations and parameters can be found in the 

State, Control, and Observation Variable Definitions section. 


CONCLUDING REMARKS 


The FORTRAN program LINEAR was developed to provide a flexible, powerful, and 
documented tool to derive linear models for aircraft stability analysis and control 
law design. This report discusses the interactive LINEAR program, describing the 
interfaces to user- supplied subroutines and input files. The output from LINEAR is 
also described. Examples of the user-supplied subroutines are presented in appen- 
dix J. A microfiche listing of the program for a VAX 11/750 (Digital Equipment 
Corporation, Maynard, Massachusetts) with the VMS operating system is included with 
this report. The program is written in standard FORTRAN 77 so that it may be porte 
to other computers without modification. 

The program LINEAR numerically determines a linear system model using nonlinear 
equations of motion and a user- supplied nonlinear aerodynamic model. LINEAR is als 
capable of extracting linearized engine effects, such as net thrust, torque, and 
gyroscopic effects and including these effects in the linear system model. The poi 
at which this linear system model is defined is determined either by completely spe 
ifying the state and control variables or by specifying an analysis point on a tra- 
jectory, selecting a trim option, and directing the program to determine the contro 
variables and remaining state variables. 


uf POOR Q 
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LINEAR has several features that make it unique among the linearization programs 
common in the aerospace industry. The most significant of these features is flexi- 
bility. By generalizing the surface definitions and making no assumptions of sym- 
metric mass distributions, the program can be applied to any aircraft in any phase 
of flight except hover. The unique trimming capability, provided by means of a user- 
supplied subroutine, allows unlimited possibilities of trimming strategies and sur- 
face scheduling, which are particularly important for asymmetrical vehicles and 
aircraft having multiple surfaces affecting a single axis. The formulation of the 
equations of motion permit the inclusion of thrust- vectoring effects. The ability 
to select, without program modification, the state, control, and observation vari- 
ables for the linear models, combined with the large number of observation quantities 
available, allows any analysis problem to be attacked with ease. 

This report documents the use of the program LINEAR, defining the equations and 
the methods used to implement the program. The trimming capabilities of LINEAR are 
discussed from both theoretical and implementation perspectives. The input and out- 
put files are described in detail. The user-supplied subroutines required for LINEAR 
are discussed and sample subroutines are presented. 


Ames Research Center 

Dry den Flight Research Facility 

National Aeronautics and Space Administration 

Edwards, California, November 3, 1987 
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APPENDIX A: CORRECTION TO AERODYNAMIC COEFFICIENTS FOR A CENTER OF GRAVITY 

NOT AT THE AERODYNAMIC REFERENCE POINT 


The point on the vehicle at which the nonlinear force and moment coefficients ai 
defined is referred to as the aerodynamic reference point. All aerodynamic effects 
are modeled at this aerodynamic reference point. Thus, when this point and the 
vehicle center of gravity are not coincident, the forces acting at the aerodynamic 
reference point effectively induce moments that act incrementally on the moments 
defined at the aerodynamic reference point by the nonlinear aerodynamic model. 

The total aerodynamic moment M acting at the vehicle center of gravity is define 
as 


M - M ar + Ar x F 


where 


«ar = [L ar M ar N ar ]T 

is the total aerodynamic moment acting at the aerodynamic reference point (denoted l 
subscript ar) of the vehicle. 


Ar = [Ax Ay Az]T 

is the displacement of the aerodynamic reference point from the vehicle center of 
gravity, and 


F = [X Y Z]T 

is the total aerodynamic force acting at the aerodynamic center where x, y, and z ai 
total forces along the x, y, and z body axes. Thus, 



L ar 

+ 

Ay 

Z - Az 

Y ‘ 

M = 

Mar 

+ 

Az 

X - Ax 

Z 


_N a r 

+ 

Ax 

Y - Ay 

X 


The total aerodynamic moment acting at the vehicle center can be expressed in 
terms of the force and moment coefficients derived from the user-supplied nonlinear 
aerodynamic modeling subroutine CCALC by defining the body axis forces in terms of 
stability axis force coefficients: 


X = qS(-C D cos a + C L sin a) 
Y = qSC Y 

Z = qS(-C D sin a - C L cos a) 
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Equations obtained by substituting these equations into the definition of the total 
I aerodynamic moment equation and applying the definitions of the total aerodynamic 
) moments, result in 

L ar = <lSbC£ ar 

! M ar = c l ScC m ar 

| N ar = qSbC naj _ 

, Expressions for total aerodynamic moment coefficients corrected to the vehicle center 
of gravity can be derived as follows: 


c £ 

C ^ar 

+ 

Ay 

b 

i 

o 

D 

sin a - 

C L 

cos a) 

Az 
“ b~ 

n 

K 

c m 

= C ">ar 

+ 

Az 

1 

n 

a 

cos a + 

C L 

sin a) 

Ax 

(-C D sin a - C L cos a) 



c 





c 



= Cn ar 

+ 

Ax 

b 

C Y ~ 

— ( C D 

cos 

a + C L 

sin 

a) 


These calculations are normally performed within LINEAR in the subroutine 
CGCALC. However, if the user elects, the calculation can be performed within the 
user-supplied aerodynamic model CCALC. 
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APPENDIX B: ENGINE TORQUE AND GYROSCOPIC EFFECTS MODEL 


Torque and gyroscopic effects represent (after thrust) the main contributions ol 
the engines to the aircraft dynamics. The torque effects arise due to thrust vec- 
tors not acting at the vehicle center of gravity. The gyroscopic effects are a con- 
sequence of the interaction of the rotating mass of the engine and the vehicle dyna- 
mics. These effects can be either major or virtually negligible, depending on the 
vehicle. 

The torque effects can be modeled by considering the thrust of an engine Fp, 

where the thrust vector is aligned with the local x axis of the engine acting at soi 
point Ar from the center of gravity of the vehicle (fig. 4). 



Figure 4. Definition of location of engine 
center of mass (cm±) relative to vehicle center 
of gravity. 


The thrust vector for the ith engine Fp^ can be defined as 

Vl f pJ t 

where F p x . , F p y .» and Fp z . are the components of thrust in the x, y, and z body 

axes, respectively. From figures 5 and 6 it can be seen that the following relatio 
ships hold: 


F p Xi = I F Pi I cos e i cos ? i 


F Pyi = I F Pi I COS e i Sin 


F P Zi “ " I F Pi Sin £ i 
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I where | Fp . | represents the magnitude of the thrust due to the ith engine, the 

angle from the thrust axis of the engine to a plane parallel to the x-y body axis 
plane, and the angle from the projection of Fp^ onto the plane parallel to the 
x-y body axis plane to the x body axis. 




Figure 5. Orientation of the engines Figure 6. Detailed definition of engine 

in the x-y and x-z body planes. location and orientation parameters. 

Denoting the point at which the thrust from the ith engine acts as Ar^, this offset 
vector can be defined as 


Afi = [Axi Ayi Az^T 


where Ax-^, Ay^, and Az^ are the x, y, and z body axis coordinates, respectively, of 
the origin of the ith thrust vector. 

The torque due to offset from the center of gravity of the ith engine At d ^ is 
then given by 


Thus, 


At. 


= Ar, x F 


Pi 



Ay 


Az 

F P yi 

Ax = 
°i 

A z 

F Pxi - 

Ax 

F P Zi 


Ax 


Ay 

F PXi 
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The total torque due to engines offset from the center of gravity of the vehicle x 0 
is given by 


n n 

T o = I At o, = I x AF Pi > 

i=1 1 i=1 

where n is the number of engines. 

For the case of vectored thrust, the equations for torque produced at the vehicl< 
center of gravity from the ith engine Ax 0 ^ are somewhat more complicated. Figure 7 

schematically represents an engine with thrust vectoring whose center of gravity is 
located at Arj_ relative to the vehicle's center of gravity. 



Figure 7 . Detailed definition of 
thrust-vectoring parameters . 


The thrust is assumed to act at -Axij in the local (engine) x axis, with the engine 

center of gravity being the origin of this local coordinate system. The thrust is 
also assumed to be vectored at angles rq and relative to the local coordinate 

axes, with nj_ being the angle from the thrust vector to the engine x-y plane and c, ^ 

the angle from the projection of the thrust vector onto the engine x-y plane to the 

local x axis. Thus, letting Fq , f' f and p' represent the x, y, and z thrust 

h x i p Yi Pz i 

components in the local engine coordinate system, repectively, these terms can be 


defined in terms of 
Ci as 

the total thrust 

for 

the 

ith 

engine Fp^ and the angles and 


F i*i ’ 1 

F 

Pi 

| cos 

^i 

cos 


F P V1 - 1 

F 

Pi 

| cos 

Hi 

sin Ci 
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sin rq 





where 


F d . = Tf' f' f' 

Pi l Pxi P yi F P Z J 


To transform this equation from the ith engine axis system to the body axis system, 
the transformation matrix 


COS £i 

cos 

5i 

-sin 

«i 

sin e.^ 

cos 

5i‘ 

cos e-^ 

cos 

«i 

cos 

?i 

sin e i 

cos 

«i 

-sin 

e i 


0 


cos 

e i 





L_ 




i 

J 


is used. 









The 

resultant 

force in 

body axis coordinates is 






F Px ± 


COS £^ COS 

-sin ^ sin 

e i 

cos 

Si 


II 

F P Yi 

= 

cos £j^ sin ^ 

cos sin 

£ i 

cos 

Si 



>i_ 


-sin e-^ 

0 

cos 

1 e i 


so that 










F Px ± = 1 

F 1 

Pi 1 

(COS Hi cos cos 

fcq COS 5 i ~ COS 

Hi 

sin 

5i 




- sin ni sin cos ) 





F Py ± = 1 

F 1 
Pi 1 

(cos n i cos ? i cos 

sin + cos 

li 

sin 

Si 


Pxi 


p Yi 


PZi 


- sxn 


hi sin cos £i ) 


F p Zi =| F Pi | (-cos rii cos ^ sin - sin cos e ± ) 

The moment arm through which the vectored thrust acts is 

" . 

A Xi - Ax Ti cos £ ^ cos 


Ari = 


Ay^ - Ax Ti cos e.^ sin 
|_Az^ + Ax t ^ sin 


and the total torque due to thrust vectoring is 


T ° = t AT o, = I * ^ Pi > 


i=1 


i=1 
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The gyroscopic effects due to the interaction of the rotating mass of the engines 
and the vehicle dynamics can be derived from the equation for the rate of change of 
angular momentum, 


T gi = h ’ei + “ X ( I e i “e i ) 

where Tg^ is the gyroscopic moment produced by the ith engine, h e ^ the rate of change 
of the angular momentum of the ith engine, I e ^ the inertia tensor for the ith engine, 
the rotational velocity of the ith engine, and co the total rotational velocity oi 
the vehicle given by 

(i) = [p q r]T 

If it is assumed that the angular momentum of the engine is constant, then 



= 0 


and the equation simplifies to 


T gi = U> x(I e± Ofe.) 

Two terms, I 6i and w ei , remain to be defined in the equation for the gyroscopic 

moment produced by the ith engine. Once again, simplifying assumptions are made. I 
is assumed that the inertia tensor of the engine contains a single nonzero entry, 

I x e . 0 O' 

0 0 0 

0 0 0 _ 

where I x is the inertia of the ith engine at the location of the ith engine, 
e i 

oriented with the local x axis coincident with the rotational velocity of the engine 
The rotational velocity of the engine has components in the x, y, and z body axes 
(p e ^, q e ^, and r ei' respectively) so that 

<%i = (Pei q ei r ei) T 

where 

p ei = | 0) ei | cos Ei cos 

q 6i = | 0) ei | cos Ei sin 


?i 

Si 
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U) ei Sin £ ± 


with 



being the total angular velocity of the ith engine. 


The engine inertia tensor must be defined in the vehicle body axis system and at 
the vehicle center of gravity. This is done in three steps. The first two steps 
involve rotating the engine inertia tensor into a coordinate system orthogonal to the 
aircraft body axis system, while the third step involves defining the engine inertia 
tensor about the vehicle center of gravity rather than about the center of gravity 
of the engine itself. First, the ith engine inertia tensor is rotated through an 
angle g-^ about the local y axis so that the new inertia tensor is oriented with its 

local x-y body axis plane parallel to the x-y body axis plane of the vehicle. The 
second step requires a rotation through an angle £i about the local z axis so that 
the local x, y, and z axes are orthogonal to the x, y, and z body axes of the 
vehicle. As determined by Gainer and Sherwood (1972) and Thelander (1965), this 
rotation is a similarity transformation that yields a new inertia tensor so 

that 


Ie i “ R ?i R £i I ei R £i R 5i 

where R^ and R e are axes transformation matrices that perform the previously 
described rotations through £ and e, respectively. These matrices are given as 


Because 


cos 0 

0 1 


-sin G-j 


sm gh 


0 


COS G 


ij 


R ?i " 


cos -sin £i 

sin cos 

0 0 


so that 


R 5i R 


i 


'cos £i cos g^ -sin cos sin e^" 

sin cos e i cos sin sin 


-sxn Si 


cos G 


i J 


-1 T 

R £i = R Si 
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and 


-1 

R Ci 


= R £ 


and the matrices are unitary. 


■1 -1 

.-.F 
'l 


T T 


R £ ,R C . = R e .R 5 . 


Rr.R e . T 

Si e l 


( R 5 i R ei) 


cos £i 
-sin 
cos £i 


i cos 


sin e^ 


sin cos 

cos 

sin sin 


-sin £^ 

0 

COS £j_ J 


Therefore , 


Ie, = I 


x ei 


cos 2 

cos 2 £i 
-cos £^ sin £j^ 


cos 2 £i cos ^i sin E, ^ -cos E-^ sin £^ cos E, 

-sin £ cos £-j_ sin E \ 


£ j_ cos 2 E,j_ 

:j_ cos sin £j_ sin 2 E,^ cos 2 £^ 

£j^ cos E,j_ -sin £^ cos £j_ sin £i 

The angular momentum of the ith engine can now be expressed as 

T 


sin^ £■■ 


h i [ h i 1 “i 2 “13 


h-: _ h 


‘J 


with 

h i 3 = Pei^xgi cos2 £ i cos2 Si] + < 2e i [ I x ei cos 2 £i cos sin 5 ± ] 

+ r ei |"-l Xe cos Ei sin e^ cos 2 f^] 

hi 2 = Pei^xgi cos2 £ i cos 5 i sin 5 i] + q e i [ I x ei sin2 Si cos2 e ] 

+ r ei|^ -I x e . sin G i cos e i sin Si] 

h i 3 = Pe i [" I x ei cos £i sin Ei cos E , ± ] + q ei [-Ix ei sin e i cos Si sin Si] 

+ r ei[lx ei sin2 £ i] 
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Thus, the gyroscopic moment induced by the ith engine x q can be expanded to 

.i 


g i 


qh;:. - rh 


13 iU i 2 


rh. 


i 1 ph i3 

LP h i 2 “ <5 h iJ 


and the total moment induced by gyroscopic interaction of the vehicle dynamics and 
the rotating engine components is 


n 



Engine torque and gyroscopic effects are modeled within the subroutine ENGINE 
using information provided by the user from the engine modeling subroutine IFENGN. 
These effects are calculated as incremental moments and are included directly in 
the equations of motion for both analysis point definition and derivation of the 
linearized system matrices. 
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APPENDIX C: OBSERVATION VARIABLE NAMES RECOGNIZED BY LINEAR 


This appendix lists all observation variable names recognized by LINEAR except 
for state and control variable names. If state variables are specified as elements 
in the observation vector, the alphanumeric descriptor must correspond to the names 
defined in appendix D of this report. When control variables are to be included in 
the observation vector, these variables must be identified exactly as they were spec- 
ified by the user. 

The input file is formatted as shown in table 1. The alphanumeric data (meas- 
urement) is left- justified in a 5A4 format. The floating point fields (PARAM) are 
used to define sensor locations not at the center of gravity. The input name spec- 
ified by the user for an observation variable serves both to identify the observation 
variable selected within the program itself as well as to identify observation 
variables on the printed output of LINEAR. 


Observation variable 

Units 

Symbol 

Alphanumeric descriptors 

Derivatives of 

state variables 

Roll acceleration 

rad/sec 2 

• 

P 

PDOT 

ROLL ACCELERATION 

Pitch acceleration 

rad/sec 2 

• 

q 

QDOT 

PITCH ACCELERATION 

Yaw acceleration 

rad/sec 2 

• 

r 

RDOT 

YAW ACCELERATION 

Velocity rate 

ft/sec 2 

• 

V 

VDOT 

VELOCITY RATE 

Angle-of-attack rate 

rad/sec 

• 

a 

ALPDOT 
ALPHA DOT 
ALPHADOT 

Angle-of-sideslip rate 

rad/sec 

• 

e 

BTADOT 
BETA DOT 
BETADOT 

Pitch attitude rate 

rad/sec 

• 

6 

THADOT 
THETA DOT 

Heading rate 

rad/sec 

• 

'i' 

PSIDOT 
PSI DOT 

Roll attitude rate 

rad/sec 

l 

PH I DOT 
PHI DOT 
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Observation variable 

Units 

Symbol 

Alphanumeric descriptors 

Derivatives 

of state 

variables 

( continued) 

Altitude rate 

ft/sec 

• 

h 

HDOT 

ALTITUDE RATE 

Velocity north 

ft/sec 

• 

X 

XDOT 

Velocity east 

ft/sec 

• 

y 

YDOT 

Accelerations 

x body axis 
acceleration 

g 

a x 

AX 

LONGITUDINAL ACCEL 
X-AXIS ACCELERATION 
X AXIS ACCELERATION 
X-BODY AXIS ACCEL 
X BODY AXIS ACCEL 

y body axis 
acceleration 

g 

a y 

AY 

Y-AXIS ACCELERATION 

Y AXIS ACCELERATION 
Y-BODY AXIS ACCEL 

Y BODY AXIS ACCEL 
LATERAL ACCELERATION 
LAT ACCEL 

LATERAL ACCEL 

z body axis 
acceleration 

g 

a z 

AZ 

Z-BODY AXIS ACCEL 
Z BODY AXIS ACCEL 

x body axis 

accelerometer at 
vehicle center of 
gravity 

g 

a nx 

ANX 

X-AXIS ACCELEROMETER 
X AXIS ACCELEROMETER 

y body axis 

accelerometer at 
vehicle center of 
gravity 

g 

a ny 

ANY 

Y-AXIS ACCELEROMETER 
Y AXIS ACCELEROMETER 

z body axis 

accelerometer at 
vehicle center of 
gravity 

g 

a nz 

ANZ 

Z-AXIS ACCELEROMETER 
Z AXIS ACCELEROMETER 


Observation variable 

Units 

Symbol 

Alphanumeric descriptors 

Accelerations 

( continued) 

normal acceleration 

g 

a n 

AN 

NORMAL ACCELERATION 
NORMAL ACCEL 
GS 
G’S 

*x body axis 

accelerometer not at 
vehicle center of 
gravity 

g 

a nx,i 

AX, I 
ANX , I 

*y body axis 

accelerometer not at 
vehicle center of 
gravity 

g 

a ny , i 

AY, I 
ANY, I 

*z body axis 

accelerometer not at 
vehicle center of 
gravity 

g 

a nz , i 

AZ,I 
ANZ, I 

*normal accelerometer 
not at vehicle 
center of gravity 

g 

a n, i 

AN, I 

load factor 

Dimension- 

less 

n 

N 

LOAD FACTOR 

Air data parameters 


Speed of sound 

ft/sec 

a 

A 

SPEED OF 

SOUND 

t Reynolds number 

Dimension- 

less 

Re 

RE 

REYNOLDS 

NUMBER 

Reynolds number 
per unit length 

ft -1 

Re' 

RE PRIME 
R/LENGTH 
R/FEET 

R/UNIT LENGTH 

Mach number 

Dimension- 

less 

M 

M 

MACH 

AMCH 



Observation variable 

Units 

Symbol 

Alphanumeric descriptors 


Air 

data parameters (continued) 

Dynamic pressure 


lb/ft 2 

q 

QBAR 

DYNAMIC PRESSURE 

Impact pressure 


lb/ ft 2 

qc 

QC 

IMPACT PRESSURE 
DIFFERENTIAL PRESSURE 

Static pressure 


lb/ft 2 

Pa 

PA 

STATIC PRESSURE 
FREESTREAM PRESSURE 

Impact/ambient 
pressure ratio 


Dimension- 

less 

‘Ic/Pa 

QC/PA 

QC/P 

Total pressure 


lb/ft 2 

Pt 

PT 

TOTAL PRESSURE 

Temperature 


deg 

absolute 

T 

TEMP 

TEMPERATURE 

FREESTREAM TEMPERATURE 

Total temperature 


deg 

absolute 

T t 

TOTAL TEMPERATURE 

Equivalent airspeed 


knot 

v e 

VEAS 

EQUIVALENT AIRSPEED 
KEAS 

Calibrated airspeed 


knot 

v c 

VCAS 

CALIBRATED AIRSPEED 
KCAS 


Flightpath-related parameters 

Flightpath angle 


rad 

Y 

GAM 

GAMMA 

FLIGHT PATH ANGLE 
GLIDE PATH ANGLE 
GLIDE SLOPE 

Flightpath 

acceleration 


g 

fpa 

FPA 

FLIGHT PATH ACCEL 

Flightpath angle 
rate 


rad/sec 

• 

Y 

GAMMA DOT 
GAMMADOT 


Observation variable 

Units 

Symbol 

Alphanumeric descriptors 

Flightpath 

-related 

parameters 

( continued) 

Vertical acceleration 

ft/sec 2 

h 

VERTICAL ACCELERATION 

H DOT DOT 

H-DOT-DOT 

HDOT-DOT 

Scaled altitude rate 

ft/sec 

h/57.3 

H-DOT/57.3 
HDOT / 57.3 

Energy-related terms 

Specific energy 

ft 

Es 

ES 

E-SUB-S 

SPECIFIC ENERGY 

Specific power 

ft/sec 

p s 

PS 

P-SUB-S 

SPECIFIC POWER 
SPECIFIC THRUST 


Force 

parameters 


Lift force 

lb 

L 

LIFT 

Drag force 

lb 

D 

DRAG 

Normal force 

lb 

N 

NORMAL FORCE 

Axial force 

lb 

A 

AXIAL FORCE 

Body axis parameters 

x body axis velocity 

ft/sec 

u 

UB 

X-BODY AXIS VELOCITY 
X BODY AXIS VELOCITY 
X-BODY AXIS VEL 
X BODY AXIS VEL 
U-BODY 
U BODY 

y body axis velocity 

ft/sec 

V 

VB 

Y-BODY AXIS VELOCITY 

Y BODY AXIS VELOCITY 
Y-BODY AXIS VEL 

Y BODY AXIS VEL 
V-BODY 

Y BODY 



Observation variable 

Units 

Symbol 

Alphanumeric descriptors 

Body axis parameters 

z body axis velocity 

ft/sec 

w 

WB 

Z-BODY AXIS VELOCITY 
Z BODY AXIS VELOCITY 
Z-BODY AXIS VEL 
Z BODY AXIS VEL 
W-BODY 
W-BODY 

Rate of change of 
velocity in 
x body axis 

ft/sec 2 

• 

u 

UBDOT 
UB DOT 

Rate of change of 
velocity in 
y body axis 

ft/sec 2 

• 

V 

VBDOT 
VB DOT 

Miscellaneous measurements not 

at vehicle center of gravity 

Rate of change of 

ft/sec 2 

• 

w 

WBDOT 

velocity in 
z body axis 



WB DOT 

*Angle of attack not 
at vehicle center 
of gravity 

rad 

a ,i 

ALPHA, I 

ALPHA INSTRUMENT 
AOA INSTRUMENT 

*Angle of sideslip not 
at vehicle center 
of gravity 

rad 

e,i 

BETA, I 

BETA INSTRUMENT 
SIDESLIP INSTRUMENT 

*Altitude instrument 
not at vehicle center 
of gravity 

ft 

h ,i 

H, I 

ALTITUDE INSTRUMENT 

*Altitude rate 

instrument not at 
vehicle center of 
gravity 

ft/sec 


HDOT , I 

Other 

miscellaneous parameters 

Vehicle total angular 
momentum 

slug-f t 2 / 
sec 2 

T 

ANGULAR MOMENTUM 
ANG MOMENTUM 


Stability axis 
rate 


roll 


rad/sec 


Ps 


STAB AXIS ROLL RATE 


Observation variable 

Units 

Symbol 

Alphanumeric descriptors 

Other miscellaneous 

parameters 

( continued) 

Stability axis pitch 
rate 

rad/sec 

qs 

STAB AXIS PITCH RATE 

Stability axis yaw 
rate 

rad/sec 

r s 

STAB AXIS YAW RATE 


*All parameters with an asterisk are for measurements at some point 
other than the vehicle center of gravity. The program LINEAR uses 
the quantities defined in the first three floating point fields as 
definitions of the location of the sensor with respect to the vehicle 
center of gravity. The three parameters define the x body, y body, 
and z body locations, in that order, of the sensor. The units of 
these offsets from the vehicle center of gravity are defined in 
units of feet. 

+ Reynolds number is defined in terms of an arbitrary unit of length 
that is input by the user. This length is input using the first 

floating point field; however, if no value is input, c is used as the 
default value. 



APPENDIX D: STATE VARIABLE NAMES RECOGNIZED BY LINEAR 


The alphanumeric descriptors specifying state variables that are recognized by 
LINEAR are listed in this appendix. In the input file, the field containing these 
descriptors uses a 5A4 format, and all characters are left justified. The input 
alphanumeric descriptor specified by the user serves both to identify the state 
variable selected by the user within the program itself as well as to identify state 
variables on the printed output of LINEAR as described in the Output Files section. 


State variable 

Units 

Symbol 

Alphanumeric descriptors 

Roll rate 

rad/sec 

P 

P 

ROLL RATE 

Pitch rate 

rad/sec 

q 

Q 

PITCH RATE 

Yaw rate 

rad/sec 

r 

R 

YAW RATE 

Velocity 

ft/sec 

V 

V 

VELOCITY 

VEL 

VTOT 

Angle of attack 

rad 

a 

ALP 

ALPHA 

ANGLE OF ATTACK 

Sideslip angle 

rad 

e 

BTA 

BETA 

SIDESLIP ANGLE 
ANGLE OF SIDESLIP 

Pitch attitude 

rad 

e 

THA 

THETA 

PITCH ATTITUDE 

Heading 

rad 

Y 

PS I 

HEADING 
HEADING ANGLE 

Roll attitude 

rad 


PHI 

ROLL ATTITUDE 
BANK ANGLE 

Altitude 

ft 

h 

H 

ALTITUDE 

Displacement north 

ft 

X 

X 

Displacement east 

ft 

y 

Y 
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APPENDIX E: ANALYSIS POINT DEFINITION IDENTIFIERS 


Analysis point definition options are selected using alphanumeric descriptors. 
These descriptors are the first record read for each analysis case. All these 
descriptors are read using a 5A4 format. The following list associates the analysis 
point definition options with their alphanumeric descriptors: 


Analysis point definition option 

Alphanumeric descriptors 

Untrimmed 

UNTRIMMED 
NO TRIM 
NONE 
NOTRIM 

Straight-and- level 

STRAIGHT AND LEVEL 
WINGS LEVEL 
LEVEL FLIGHT 

Pushover-pul lup 

PUSHOVER AND PULLUP 
PULLUP 

PUSH-OVER/PULL-UP 
PUSH-OVER / PULL-UP 
PUSHOVER / PULLUP 
PUSHOVER PULLUP 
PUSH OVER 
PULL UP 
PUSHOVER 
PUSHPULL 

Level turn 

LEVEL TURN 


WINDUP TURN 

Thrust stabilized turn 

THRUST STABILIZED TURN 
THRUST LIMITED TURN 
FIXED THROTTLE TURN 
FIXED THRUST TURN 

Beta 

BETA 


SIDESLIP 

Specific power 

SPECIFIC POWER 
PS 

P-SUB-S 


Each of these analysis point definitions except the untrimmed, beta, and specific 
power options has two suboptions associated with it. The suboptions are requested 
using alphanumeric descriptors read using an A4 format. These suboptions are defined 
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in the Analysis Point Definition section of this report. The following list defines 
these suboptions and the alphanumeric descriptors associated with each: 


Analysis point definition 

Alphanumeric descriptors 

suboptions 


Straight-and-level : 


Alpha-trim 

ALP 

ALPH 

ALPHA 

Mach-trim 

MACH 


AMCH 

Pushover-pullup : 


Alpha-trim 

ALP 


ALPH 


ALPHA 

Load-factor-trim 

LOAD 


GS 


G’S 


AN 

Level turn: 


Alpha-trim 

ALP 

ALPH 

ALPHA 

Load-f actor- trim 

LOAD 


GS 


G'S 


AN 

Thrust stabilized turn: 


Alpha-trim 

ALP 

ALPH 

ALPHA 

Load-f actor- trim 

LOAD 


GS 


G'S 


AN 
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APPENDIX F: EXAMPLE INPUT FILE 


The following listing is an example of an input file to LINEAR. This file was 
used with the example subroutines listed in appendix J to generate the analysis and 
line printer files shown in appendixes H and I, respectively. The formats specified 
in table 1 would be the same whether each type of input has its own file or all of 
the data is in one file. 

USER'S GUIDE 

6 . 080000E+02 4.280000E+01 1.595000E+01 4.500000E+04 
2 . 870000E+04 1.651000E+05 1 .879000E+05-5.200000E+02 0. 0. 

0. 0. 0. CCALC WILL CALCULATE CG CORRECTIONS 

-1 .000000E+01 4 . 000000E+0 1 

4 STAN RADI 

ALPHA 

Q 

THETA 

VEL 

3 

ELEVATOR 5 

THROTTLE 1 2 

SPEED BRAKE 10 

2 STAN 
AN 
AY 

-2 . 900E+00 5.430E+00-4 . 000E+00 4.000E+00-3.250E+00 3.250E+00-1 .000E+00 1.000E+00 
0 ADDITIONAL SURFACES 


WINDUP TURN 

ALPHA 

H 

20000.0 

MACH 

0.90 

AN 

3.00 

BETA 

O 

o 

NEXT 

LEVEL FLIGHT 

ALPHA 

H 

20000.0 

MACH 

0.9 

GAMMA 

10.0 

END 



This input file is for a case called USER'S GUIDE (record 1). The project title 
is input from the terminal as described in the Data Input section. Record 2 speci- 
fies the mass and geometric properties of the vehicle as 

S = 608 ft 2 

b = 42.8 ft 
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c = 15.95 ft 
w = 45,000 lbs 

Record 3 defines the moments and products of intertia of the vehicle as 


X 

= 28,700 

slug-ft 2 

y 

= 165,100 

slug-ft 2 

z 

= 187,900 

slug-ft 2 

xz 

= -520 

slug-ft 2 

xy 

= 0.0 

slug-f t 2 

yz 

= 0.0 

slug-f t 2 


Record 4 defines the location of the aerodynamic reference point to be coincident 
with the vehicle center of gravity of the nonlinear aerodynamic model by setting 


Ax = 0.0 ft 


Ay = 0.0 ft 
Az = 0.0 ft 

Record 4 also specifies that LINEAR should not use its internal model to make 
corrections for the offset in the vehicle center of gravity from the aerodynamic 
reference point because the aerodynamic model includes such corrections. 

Record 5 defines the angle-of-attack range of the aerodynamic model. 

Record 6 specifies that there will be four state variables in the linear 
model, that the formulation of the state equation will be in the standard form 

(x = A x + B u), and that the nondimensional stability and control derivatives with 

respect to angle of attack and angle of sideslip should be scaled in radians -1 . The 
next four records define the output formulation of the state vector to be 


x 


“a” 

q 

0 

_v_ 


Record 11 specifies that the linear model will have 
control vector. The following three records (12 to 14) 


three parameters in the 
specify that 


u 


elevator 
throttle 
speed brake 
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and that elevator, throttle, and speed brake are located in DC(5), DC(12), and DC(10) 
of the /CONTROL/ common block. 


Record 15 specifies that two observation variables will be used, and that the 
observation equation will be in the standard form, 

y = Hx + Fu 

The next two records (16 and 17) define the elements of the output vector to be 



Record 18 specifies the ranges for the trim parameters DES, DAS, DRS , and THRSTX , 
used to trim the longitudinal, lateral, and directional axes, and thrust, respec- 
tively. The ranges for these parameters are defined by record 20 to be 


-2.9 < DES <5.43 
-4.0 < DAS <4.0 
-3.25 < DRS < 3.25 
-1.0 < THRSTX < 1.0 


The first three parameters essentially represent stick and rudder positions and are 
so specified because of the implementation of the subroutine UCNTRL (discussed in 
app. J). The thrust trim parameter is specified in this manner because of the imple- 
mentation of UCNTRL. Speed brake is scheduled when THRSTX < 0 and thrust is com- 
manded when THRSTX > 0. 

Record 19 specifies that no additional control surfaces are to be set. 

The next seven records (20 to 26) define an analysis point option. These records 
request a level turn trim option at 


h = 20,000 ft 
M = 0.9 
a n = 3.0 g 


and 


0 = 0 . 0 ° 

The second record of this set (record 21) indicates which level turn suboption is 
requested. The alphanumeric descriptor ALPHA indicates that angle of attack is to be 
varied until the specified 3.0-g turn is achieved. The final record of this analysis 
point option definition set contains the keyword NEXT to indicate both an end to the 
current analysis point option definition and that another analysis point option defi- 
nition follows. 
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The final six records {records 27 to 32) define a straight-and-level analysis 
point option at 


h = 20,000 ft 
M = 0.9 


and 


y = 10.0° 

The second record of this set (record 28) identifies the Alpha Trim suboption in 
which angle of attack is varied until trim is achieved at the specified condition. 
The final record of this set contains the keyword END to indicate the termination of 
the current analysis point definition as well as the termination of input cases. 
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APPENDIX G: EXAMPLE LINEARIZED STABILITY AND CONTROL DERIVATIVE FILE 


STABILITY AND CONTROL DERIVATIVES FOR CASE #1 
LINEARIZER TEST AND DEMONSTRATION CASES 


FOR THE PROJECT: 

F-15 SIMULATION 

REFERENCE ALTITUDE = 2.0000D+04 

REFERENCE MACH = 9.0000D-01 

ANGLES ARE IN RADIANS FOR OUTPUT 

WARNING!!!! DERIVATIVES WITH RESPECT TO MACH WILL BE USED. 
VELOCITY DERIVATIVES WILL BE COMPUTED FROM MACH DERIVATIVES. 


THERE ARE 3 CONTROL SURFACES 


ROLL MOMENT DERIVATIVES 


CLO 

ROLL RATE 
PITCH RATE 
YAW RATE 
VELOCITY 
MACH NUMBER 
ALPHA 
BETA 
ALTITUDE 
ALPHA DOT 
BETA DOT 


(RAD/SEC) 
(RAD/SEC) 
( RAD/SEC ) 
(FT/SEC) 

(RAD) 
(RAD) 
(FT) 
( RAD/SEC ) 
(RAD/SEC) 


ELEVATOR 
THROTTLE 
SPEED BRAKE 


5 

12 

10 


PITCH MOMENT DERIVATIVES 


CMO 

ROLL RATE (RAD/SEC) 
PITCH RATE (RAD/SEC) 
YAW RATE ( RAD/SEC ) 

VELOCITY (FT/SEC) 
MACH NUMBER 
ALPHA 
BETA 
ALTITUDE 
ALPHA DOT 
BETA DOT 

ELEVATOR 
THROTTLE 
SPEED BRAKE 


(RAD) 
(RAD) 
(FT) 
(RAD/SEC) 
(RAD /SEC) 


5 

12 

10 


-4.02966D-05 
-2 . OQOOOD-O 1 
. 00000D+00 
1 .50990D-01 
-1.26696D-07 
-1 .31375D-04 
.OOOOOD+OO 
-1 .33450D-0 1 
.OOOOOD+OO 
•OOOOOD+OO 
.OOOOOD+OO 

.00000D+00 

.OOOOOD+OO 

.OOOOOD+OO 


4.22040D-02 
.OOOOOD+OO 
3.89527D+00 
.OOOOOD+OO 
-3 . 28490D-06 
-3 . 40620D-03 
-1.68819D-01 
.OOOOOD+OO 
.00000D+00 
-1 . 18870D+0 1 
.00000D+00 

-6.9528 1D-0 1 
.OOOOOD+OO 
-4. 17500D-0 1 
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YAW MOMENT DERIVATIVES 


CNO 

ROLL RATE ( RAD/SEC ) 
PITCH RATE (RAD/SEC) 
YAW RATE ( RAD/SEC ) 
VELOCITY (FT/SEC) 
MACH NUMBER 
ALPHA 
BETA 
ALTITUDE 
ALPHA DOT 
BETA DOT 

ELEVATOR 
THROTTLE 
SPEED BRAKE 

DRAG FORCE DERIVATIVES 


(RAD) 

(RAD) 

(FT) 

(RAD/SEC) 

(RAD/SEC) 


CDO 

ROLL RATE ( RAD/SEC ) 
PITCH RATE (RAD /SEC) 
YAW RATE (RAD/SEC) 
VELOCITY (FT/SEC) 
MACH NUMBER 
ALPHA 
BETA 

ALTITUDE 
ALPHA DOT 
BETA DOT 

ELEVATOR 
THROTTLE 
SPEED BRAKE 

LIFT FORCE DERIVATIVES 


(RAD) 
(RAD) 
(FT) 
(RAD/SEC) 
( RAD/SEC ) 


CLFTO 


ROLL RATE (RAD/SEC) 
PITCH RATE (RAD/SEC) 
YAW RATE ( RAD/SEC ) 
VELOCITY (FT/SEC) 
MACH NUMBER 


ALPHA 
BETA 
ALTITUDE 
ALPHA DOT 
BETA DOT 


(RAD) 
(RAD) 
(FT) 
( RAD/SEC ) 
( RAD/SEC ) 


ELEVATOR 
THROTTLE 
SPEED BRAKE 


2.25747D-04 
-3.37217D-02 
•OOOOOD+OO 
-4.0471 OD-O 1 
3.21 400D-07 
3 . 33268D-04 
. 00000D+00 
1 .29960D-01 
. 00000D+00 
. 00000D+00 
. 00000D+00 

5 . 00000D+00 

12 . 00000D+00 

10 .00000D+00 


1 .08760D-02 
. 00000D+00 
. OOOOOD+OO 
. 00000D+00 
.00000D+00 
.00000D+00 
3 . 72570D-0 1 
.00000D+00 
.00000D+00 
•00000D+00 
.00000D+00 

5 4 . 383 18D-02 

12 .00000D+00 

10 6 . 49355D-02 


1 .57360D-01 
.00000D+00 
-1.7231 5D+0 1 
.00000D+00 
1 .45286D-05 
1 .50651D-02 
4.87061D+00 
.00000D+00 
.OOOOOD+OO 
1 .72315D+01 
.00000D+00 

5 5 . 72950D-0 1 

12 .00000D+00 

10 3 . 749 13D-02 



SIDE FORCE DERIVATIVES 


CYO 



5 . 32725D-04 

ROLL RATE 

{ RAD/SEC ) 


.00000D+00 

PITCH RATE 

(RAD/SEC) 


. 00000D+00 

YAW RATE 

(RAD/SEC) 


• 00000D+00 

VELOCITY 

(FT/SEC) 


• 00000D+00 

MACH NUMBER 



. 00000D+00 

ALPHA 

(RAD) 


.00000D+00 

BETA 

(RAD) 


-9.74030D-01 

ALTITUDE 

(FT) 


. 00000D+00 

ALPHA DOT 

(RAD/SEC) 


. 00000D+00 

BETA DOT 

(RAD/SEC) 


. 00000D+00 

ELEVATOR 


5 

. OOOOOD+OO 

THROTTLE 


12 

•OOOOOD+OO 

SPEED BRAKE 


10 

•00000D+00 
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STABILITY AND CONTROL DERIVATIVES FOR CASE #2 
LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

F-15 SIMULATION 

REFERENCE ALTITUDE = 2.0000D+04 

REFERENCE MACH = 9.0000D-01 

ANGLES ARE IN RADIANS FOR OUTPUT 

WARNING!!!! DERIVATIVES WITH RESPECT TO MACH WILL BE USED. 
VELOCITY DERIVATIVES WILL BE COMPUTED FROM MACH DERIVATIVES. 

THERE ARE 3 CONTROL SURFACES 


ROLL MOMENT DERIVATIVES 


CLO 



1 .25377D-16 

ROLL RATE 

( RAD/SEC ) 


-2 . OOOQOD-O 1 

PITCH RATE 

(RAD/SEC) 


. OOOOOD+OO 

YAW RATE 

(RAD/SEC) 


1 .50990D-01 

VELOCITY 

(FT/SEC) 


. 00000D+00 

MACH NUMBER 


.OOOOOD+OO 

ALPHA 

(RAD) 


. 00000D+0Q 

BETA 

(RAD) 


-1 .33450D-01 

ALTITUDE 

(FT) 


.OOOOOD+OO 

ALPHA DOT 

(RAD/SEC) 


.OOOOOD+OO 

BETA DOT 

( RAD/SEC ) 


.00000D+00 

ELEVATOR 


5 

.OOOOOD+OO 

THROTTLE 


12 

.OOOOOD+OO 

SPEED BRAKE 


10 

.OOOOOD+OO 

PITCH MOMENT DERIVATIVES 



CMO 



4 . 22040D-02 

ROLL RATE 

(RAD/SEC) 


.OOOOOD+OO 

PITCH RATE 

( RAD/SEC ) 


3 . 89530D-00 

YAW RATE 

(RAD/SEC) 


.OOOOOD+OO 

VELOCITY 

( FT/SEC ) 


-7 . 79749D- 1 0 

MACH NUMBER 



-8.08541D-07 

ALPHA 

(RAD) 


-1.68819D-01 

BETA 

(RAD) 


.OOOOOD+OO 

ALTITUDE 

(FT) 


.OOOOOD+OO 

ALPHA DOT 

(RAD/SEC) 


-1 . 18870D+0 1 

BETA DOT 

(RAD/SEC) 


•OOOOOD+OO 

ELEVATOR 


5 

-6 . 95279D-0 1 

THROTTLE 


12 

.00000D+00 

SPEED BRAKE 


10 

-4 . 1 7500D-0 1 
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YAW MOMENT DERIVATIVES 


CNO 



1 .22535D-16 

ROLL RATE 

( RAD/SEC ) 


-3.3721 7D-02 

PITCH RATE 

(RAD/SEC) 


.00000D+00 

YAW RATE 

{ RAD/SEC) 


-4.04710D-01 

VELOCITY 

( FT/SEC ) 


. 00000D+00 

MACH NUMBER 


. Q0000D+00 

ALPHA 

(RAD) 


•OOOOOD+OO 

BETA 

(RAD) 


1 .29960D-0 1 

ALTITUDE 

(FT) 


.OOOOOD+OO 

ALPHA DOT 

(RAD/SEC) 


•00000D+00 

BETA DOT 

(RAD/SEC) 


•OOOOOD+OO 

ELEVATOR 


5 

.00000D+00 

THROTTLE 


12 

.00000D+00 

SPEED BRAKE 


10 

.00000D+00 

DRAG FORCE 

DERIVATIVES 



CDO 



1.08760D-02 

ROLL RATE 

(RAD/SEC) 


•00000D+00 

PITCH RATE 

(RAD/SEC) 


.00000D+00 

YAW RATE 

(RAD/SEC) 


•OOOOOD+OO 

VELOCITY 

(FT/SEC) 


•OOOOOD+OO 

MACH NUMBER 



•OOOOOD+OO 

ALPHA 

(RAD) 


3.72570D-01 

BETA 

(RAD) 


.00000D+00 

ALTITUDE 

(FT) 


•OOOOOD+OO 

ALPHA DOT 

(RAD/SEC) 


.00000D+00 

BETA DOT 

(RAD/SEC) 


.00000D+00 

ELEVATOR 


5 

4 . 383 1 3D-02 

THROTTLE 


12 

•OOOOOD+OO 

SPEED BRAKE 


10 

6 . 49346D-02 

LIFT FORCE 

DERIVATIVES 



CLFTO 



1 . 57360D-0 1 

ROLL RATE 

(RAD/SEC) 


•OOOOOD+OO 

PITCH RATE 

(RAD/SEC) 


-1 .72320D+01 

YAW RATE 

(RAD/SEC) 


.00000D+00 

VELOCITY 

(FT/SEC) 


3 . 59263D-09 

MACH NUMBER 



3.72529D-06 

ALPHA 

(RAD) 


4. 87061 D+ 00 

BETA 

(RAD) 


.00000D+00 

ALTITUDE 

(FT) 


.00000D+00 

ALPHA DOT 

(RAD/SEC) 


1 .72320D+01 

BETA DOT 

(RAD/SEC) 


•OOOOOD+OO 

ELEVATOR 


5 

5 . 7296 1D-0 1 

THROTTLE 


12 

.00000D+00 

SPEED BRAKE 


10 

3.7491 3D-02 


88 



SIDE FORCE DERIVATIVES 


CYO 

ROLL RATE ( RAD/SEC ) 
PITCH RATE (RAD/SEC) 
YAW RATE ( RAD/SEC ) 

VELOCITY (FT/SEC) 
MACH NUMBER 
ALPHA 
BETA 

ALTITUDE 
ALPHA DOT 
BETA DOT 

ELEVATOR 
THROTTLE 
SPEED BRAKE 


4. 15737D-16 
.OOOOOD+OO 
.OOOOOD+OO 
. 00000D+00 
.OOOOOD+OO 
.00000D+00 
.OOOOOD+OO 
-9 . 74030D-0 1 
•OOOOOD+OO 
.OOOOOD+OO 
.OOOOOD+OO 

5 .OOOOOD+OO 

12 .OOOOOD+OO 

10 .OOOOOD+OO 


(RAD) 
(RAD) 
(FT) 
( RAD/SEC ) 
(RAD/SEC) 
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APPENDIX H: EXAMPLE OUTPUT ANALYSIS FILE 


The following is an example analysis file. This file was produced using the 
example input file listed in appendix F and the example user-supplied subroutines 
listed in appendix J. 


LINEARIZER TEST AND DEMONSTRATION CASES 
USER'S GUIDE 


TEST CASE # 




X - DIMENSION =4 U - DIMENSION =3 Y - DIMENSION = 

STATE EQUATION FORMULATION: STANDARD 

OBSERVATION EQUATION FORMULATION: STANDARD 


1 

2 


STATE VARIABLES 
ALPHA 

Q 

THETA 

VEL 


0 . 465695D-0 1 
0 . 92 1 683D-0 1 
0. 159885D-01 
0 . 933232D+03 


RADIANS 

RADIANS/SECOND 

RADIANS 

FEET/SECOND 


CONTROL VARIABLES 


ELEVATOR 
THROTTLE 
SPEED BRAKE 


= 0 . 538044D-0 1 
= 0.2141 05D+00 
= 0 . 000000D+00 


DYNAMIC INTERACTION VARIABLES 


X-BODY AXIS FORCE 
Y-BODY AXIS FORCE 
Z-BODY AXIS FORCE 
PITCHING MOMENT 
ROLLING MOMENT 
YAWING MOMENT 


0. 134742D+06 
0. 102657D+05 
0 . 9381 55D-03 
0.680525D+02 
0 . 000000D+00 
- . 2 1 200 7D+02 


POUNDS 

POUNDS 

POUNDS 

FOOT-POUNDS 

FOOT-POUNDS 

FOOT-POUNDS 


OBSERVATION VARIABLES 


AN 

AY 


0 . 300 163D+0 1 GS 

0.941 435D+00 GS 
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A-MATRIX FOR: DX / DT 


A * X + B * U + D * V 


-0. 121436D+01 0 . 1 00000D+0 1 0 . 1 36756D-02-0 . 1 2 1605D-03 
-0. 147423D+01-0.221451D+01-0.450462D-02 0.294019D-03 
0 . 000000D+00 0.331812D+00 0.000000D+00 O.OOOOOOD+OO 
-0.790853D+02 0 . OOOOOOD+OO-O . 320822D+02-0 . 157297D-0 1 

B-MATRIX FOR: DX/DT = A*X + B * U + D * V 

0. 141961D+00-0. 164948D-02-0 . 928933D-02 
0 .220778D+02 0 . 543324D-02-0 . 135074D+02 
0.000000D+00 0.000000D+00 0.000000D+00 
0. 105186D+02 0 . 34281 7D+02-0 . 1 55832D+02 


D-MATRIX FOR: DX/DT = A*X + B*U + D*V 


0 . 343642D-07 0.000000D+00 0.737378D-06 O.OOOOOOD+OO 0.000000D+00 


0.11 3 192D-06 0. 000000D+00-0. 242885D-05 0.605694D-05 O.OOOOOOD+OO 
O.OOOOOOD+OO O.OOOOOOD+OO O.OOOOOOD+OO O.OOOOOOD+OO O.OOOOOOD+OO 
0 . 7 14203D-03 0.398492D-06 0.332842D-04 O.OOOOOOD+OO O.OOOOOOD+OO 


O.OOOOOOD+OO 

O.OOOOOOD+OO 

O.OOOOOOD+OO 

O.OOOOOOD+OO 


H-MATRIX FOR: Y = H*X + F*U+E*V 

0 . 35 1 752D+02 0.277556D-16 0.150046D-02 0.640771D-02 
O.OOOOOOD+OO 0. 000000D+00-0. 150534D-01 0.446789D-10 


F-MATRIX FOR: Y = H*X + F*U + E*V 

0.41 2845D+0 1 -0 . 180978D-02 0.291699D+00 
O.OOOOOOD+OO O.OOOOOOD+OO O.OOOOOOD+OO 


E-MATRIX FOR: Y=H*X+F*U+E*V 


-0.377037D-07 0. 000000D+00-0. 2 14 132D-04 
O.OOOOOOD+OO 0 . 222222D-04 O.OOOOOOD+OO 
O.OOOOOOD+OO 
O.OOOOOOD+OO 


O.OOOOOOD+OO O.OOOOOOD+OO 
O.OOOOOOD+OO O.OOOOOOD+OO 


TEST CASE *************** ******************************** * 

X - DIMENSION =4 U - DIMENSION =3 Y - DIMENSION = 

STATE EQUATION FORMULATION: STANDARD 

OBSERVATION EQUATION FORMULATION: STANDARD 


2 

2 


STATE VARIABLES 
ALPHA 

Q 

THETA 

VEL 


= - . 1 26650D-0 1 
= 0 . 000000D+00 
= 0. 161868D+00 
= 0.933232D+03 


RADIANS 

RADIANS/SECOND 

RADIANS 

FEET/SECOND 


CONTROL VARIABLES 


ELEVATOR 
THROTTLE 
SPEED BRAKE 


= 0.637734D-01 
= 0 . 225092D+00 
= 0 . 000000D+00 


DYNAMIC INTERACTION VARIABLES 


X-BODY AXIS FORCE 
Y-BODY AXIS FORCE 
Z-BODY AXIS FORCE 
PITCHING MOMENT 
ROLLING MOMENT 
YAWING MOMENT 


0.443769D+05 
0.300494D+04 
0.445795D-09 
-.528481D-1 1 
0 . 000000D+00 
0. 114239D-10 


POUNDS 

POUNDS 

POUNDS 

FOOT-POUNDS 

FOOT-POUNDS 

FOOT-POUNDS 


OBSERVATION VARIABLES 

AN = 0.985228D+00 GS 

AY = 0.990655D-14 GS 


A-MATRIX FOR: DX/DT = A*X + B*U+D*V 

0. 120900D+01 0. 100000 D+01-0.575730D-02-0. 70 1975D-04 
0 . 1 49 189D+0 1-0.221 45 1D+0 1 0.189640D-01 0.231368D-03 
0 . OOOOOOD+OO 0 . 100000D+0 1 0.000000D+00 O.OOOOOOD+OO 
0 . 576868D+02 0 . OOOOOOD+OO-O . 3 1 625 1D+02-0 .460435D-02 


B-MATRIX FOR: DX/DT = A*X + B * U + D * V 

0. 141961D+00 0 • 4487 42D-0 3-0 . 9 289 3 2D- 02 
0.220778D+02-0. 14781 2D-02-0. 135074D+02 
0.000000D+00 O.OOOOOOD+OO 0.000000D+00 
0. 105186D+02 0 . 343 162D+02-0 . 155832D+02 
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\ D-MATRIX FOR: DX/DT+A*X+B*U+D*V 

0.934880D-08 0.000000D+00 0.738119D-06 0.000000D+00 0.000000D+00 
0 . 30794 ID-07 0 . OOOOOOD+OO-O . 243 1 29D-05 0.605694D-05 0.000000D+00 
0 . 000000D+00 0 . 000000D+00 0.000000D+00 0.000000D+00 0.000000D+00 
0 . 714920D-03 0 . OOOOOOD+OO-O . 905497D-05 0.000000D+00 0.000000D+00 
0.000000D+00 
0.000000D+00 
0 . 000000D+00 
0 . 000000D+00 


H-MATRIX FOR: Y = H*X + F*U 

0.350424D+02 0 . 277556D-16-0 . 6323 14D-02 0.203434D-02 
0 . 000000D+00 0.000000D+00 O.OOOOOOD+OO 0.212306D-16 

F-MATRIX FOR: Y = H*X + F*U 

0.41 1323D+01 0.492845D-03 0.263288D+00 
0.000000D+00 O.OOOOOOD+OO 0.000000D+00 


E-MATRIX FOR: Y+H*X+F*U+E*V 

0 . 1 02676D-07 0 . OOOOOOD+OO-O . 2 141 16D-04 O.OOOOOOD+OO O.OOOOOOD+OO 
O.OOOOOOD+OO 0 . 222222D-04 O.OOOOOOD+OO O.OOOOOOD+OO O.OOOOOOD+OO 
O.OOOOOOD+OO 
O.OOOOOOD+OO 




APPENDIX I: EXAMPLE PRINTER OUTPUT FILE 


The following listings are printer output files generated by LINEAR using the 
example input file in appendix F and the example user-supplied subroutines in appen- 
dix J. 


GEOMETRY AND MASS DATA FOR: 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 


WING 

AREA 

608.000 

(FT**2 ) 

WING 

SPAN 

42.800 

(FT) 

MEAN 

CHORD 

15.950 

(FT) 

VEHICLE WEIGHT 

45000.000 

(LB) 

IX 


28700.000 

( SLUG— FT**2 ) 

IY 


165100.000 

(SLUG-FT**2 ) 

IZ 


187900.000 

(SLUG— FT* *2) 

IXZ 


-520.000 

( SLUG-FT**2 ) 

IXY 


0.000 

(SLUG-FT **2) 

IYZ 


0.000 

( SLUG— FT* *2 ) 


VECTOR DEFINING REFERENCE POINT OF AERODYNAMIC MODEL 
WITH RESPECT TO VEHICLE CENTER OF GRAVITY: 

DELTA X 0.000 (FT) 

DELTA Y 0.000 (FT) 

DELTA Z 0.000 (FT) 


FORCE AND MOMENT COEFFICIENT CORRECTIONS FROM THE REFERENCE POINT OF 
THE BASIC AERODYNAMIC MODEL TO THE VEHICLE CENTER OF GRAVITY ARE 
INTERNALLY CALCULATED IN LINEAR. 

MINIMUM ANGLE OF ATTACK -10.000 (DEG) 

MAXIMUM ANGLE OF ATTACK 40.000 (DEG) 


PARAMETERS USED IN THE STATE VECTOR FOR: 
LINEARIZER TEST AND DEMONSTRATION CASES 


94 





FOR THE PROJECT: 

USER'S GUIDE 

ALPHA 

Q 

THETA 

VEL 

THE STANDARD FORMULATION OF THE STATE EQUATION 
HAS BEEN SELECTED. THE FORM OF THE EQUATION IS: 

DX/DT =A*X+B*U 

SURFACES TO BE USED FOR THE CONTROL VECTOR FOR: 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER’S GUIDE 

LOCATION IN /CONTROL/ 

ELEVATOR 5 

THROTTLE 12 

SPEED BRAKE 1 0 


PARAMETERS USED IN THE OBSERVATION VECTOR FOR: 
LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 

AN 

AY 


THE STANDARD FORMULATION OF THE OBSERVATION EQUATION 
HAS BEEN SELECTED. THE FORM OF THE EQUATION IS: 

Y = H*X + F*U 


LIMITS FOR TRIM OUTPUT PARAMETERS: 

MINIMUM MAXIMUM 

PITCH AXIS PARAMETER -2.900 5.430 

ROLL AXIS PARAMETER -4.000 4.000 

YAW AXIS PARAMETER -3.250 3.250 

THRUST PARAMETER -1.000 1.000 


NO ADDITIONAL SURFACES TO BE SET WERE DEFINED 


TRIM CONDITIONS FOR CASE # 1 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 

LEVEL TURN WHILE VARYING ALPHA 
TRIM ACHIEVED 


COEFFICIENT OF LIFT 


= 

0.40144 

COEFFICIENT OF DRAG 


= 

0.03058 

LIFT 

(LBS) 

= 

134741.68807 

DRAG 

(LBS) 

= 

10265.70661 

ALTITUDE 

(FT) 

= 

20000.00000 

MACH 


= 

0.90000 

VELOCITY 

(FT/SEC) 

= 

933.23196 

EQUIVALENT AIRSPEED 

(KTS ) 

= 

403.42303 

SPEED OF SOUND 

(FT/ SEC) 

= 

1036.92440 

GRAVITATIONAL ACCEL 

( FT/ SEC**2 ) 

= 

32.11294 

NORMAL ACCELERATION 

(G-S) 

= 

3.00163 

LOAD FACTOR 


= 

2.99995 

DYNAMIC PRESSURE 

(LBS/FT**2 ) 

= 

552.05302 

DENSITY 

( SLUG/FT**3 ) 

= 

0.00126774 

WEIGHT (^ALTITUDE) 

(LBS) 

= 

44914.60434 

BETA 

(DEG) 

= 

0.03193 

ALPHA 

(DEG) 

= 

2.66824 

PHI 

(DEG) 

= 

70.62122 

THETA 

(DEG) 

= 

0.91607 

ALTITUDE RATE 

(FT/SEC) 

= 

0.00000 

GAMMA 

(DEG) 

= 

0.00000 

ROLL RATE 

(DEG/SEC) 

= 

-0.08951 

PITCH RATE 

( DEG/SEC ) 

= 

5.28086 

YAW RATE 

(DEG/SEC) 

= 

1.85749 

THRUST 

(LBS) 

= 

10277.03515 

SUM OF THE SQUARES 


= 

0.00000 

TRIM PARAMETERS 

TRIM PITCH AXIS PARAMETER 

2= 

-0.66958 

TRIM ROLL AXIS PARAMETER 

= 

-0.01526 

TRIM YAW AXIS PARAMETER 

= 

-0.02125 

TRIM THRUST PARAMETER 


= 

0.21410 

CONTROL VARIABLES 

ELEVATOR 


= 

0.05380 

THROTTLE 


= 

0.21410 

SPEED BRAKE 

(DEG) 

= 

0.00000 

OBSERVATION VARIABLES 

AN 



3.00163339 

AY 


= 

.94136286 
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NON-DIMENSIONAL STABILITY AND CONTROL DERIVATIVES FOR CASE 
LINEARIZER TEST AND DEMONSTRATION CASES 
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MATRIX A USING THE FORMULATION OF THE STATE EQUATION: 

DX/DT = A * X + B * U + D * V 
FOR CASE # 1 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 

ALPHA Q THETA VEL 

TIME 

DERIVATIVES 
WRT : 

ALPHA -0 . 1 2 1436D+0 1 0.100000D+01 0.136756D-02 -0 • 12 1 605D-03 

Q -0. 147423D+01 -0 . 22 145 1D+0 1 -0 .450462D-02 0.294019D-03 

THETA 0 . OOOOOOD+OO 0.331812D+00 0.000000D+00 0.000000D+00 

VEL -0.790853D+02 0. OOOOOOD+OO -0 .320822D+02 -0 . 1 57297D-0 1 


MATRIX B USING THE FORMULATION OF THE STATE EQUATION: 
DX/DT = A * X + B * U + D * V 
FOR CASE # 1 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 


ELEVATOR THROTTLE SPEED BRAKE 

TIME 

DERIVATIVES 

WRT: 

ALPHA -0. 141961D+00 -0. 164948D-02 -0 . 928933D-02 

Q -0 . 220778D+02 0.543324D-02 -0 . 135074D+02 

THETA 0. 000000D+00 0. OOOOOOD+OO 0. 000000D+00 

VEL -0. 105186D+02 0.342817D+02 -0. 155832D+02 
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MATRIX D USING THE FORMULATION OF THE STATE EQUATION 
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MATRIX H USING THE FORMULATION OF THE OBSERVATION EQUATION: 


Y = H * X + F * U + E * V 

FOR CASE # 1 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER’S GUIDE 

ALPHA Q THETA VEL 

AN 0 . 35 1 752D+02 0.277556D-16 0.150046D-02 0.640771D-02 

AY 0 . 000000D+00 O.OOOOOOD+OO -0 . 1 50534D-0 1 0.446789D-10 


MATRIX F USING THE FORMULATION OF THE OBSERVATION EQUATION: 
Y = H * X + F * U + E * V 

FOR CASE # 1 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 

ELEVATOR THROTTLE SPEED BRAKE 

AN 0 . 4 12845D+0 1 -0. 180978D-02 0.291699D+00 

AY O.OOOOOOD+OO O.OOOOOOD+OO 0.000000D+00 
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MATRIX E USING THE FORMULATION OF THE OBSERVATION EQUATION 
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TRIM CONDITIONS FOR CASE # 2 

LINEARIZER TEST AND DEMONSTRATION CASES 


FOR THE PROJECT: 

USER'S GUIDE 

STRAIGHT AND LEVEL TRIM WHILE VARYING ALPHA 
TRIM ACHIEVED 


COEFFICIENT OF LIFT 
COEFFICIENT OF DRAG 
LIFT 
DRAG 

ALTITUDE 

MACH 

VELOCITY 

EQUIVALENT AIRSPEED 
SPEED OF SOUND 
GRAVITATIONAL ACCEL 
NORMAL ACCELERATION 
LOAD FACTOR 
DYNAMIC PRESSURE 
DENSITY 

WEIGHT (^ALTITUDE) 

BETA 

ALPHA 

PHI 

THETA 

ALTITUDE RATE 
GAMMA 
ROLL RATE 
PITCH RATE 
YAW RATE 
THRUST 

SUM OF THE SQUARES 



= 

0 


= 

0 

(LBS) 

= 

44376 

(LBS) 

= 

3004 

(FT) 

= 

20000 


= 

0 

(FT/SEC) 

= 

933 

(KTS ) 

= 

403 

( FT/SEC ) 

= 

1036 

(FT/SEC* *2) 

= 

32 

(G-S) 

= 

0 


= 

0 

(LB/FT* * 2 ) 

= 

552 

( SLUG/FT**3 ) 

= 

0 

(LBS) 

= 

44914 

(DEG) 

= 

0 

(DEG) 

= 

-0 

(DEG) 

= 

0 

(DEG) 

= 

9 

(FT/SEC) 

= 

162 

(DEG) 

= 

10 

(DEG/SEC) 

= 

0 

(DEG/SEC) 

= 

0 

(DEG/SEC) 

= 

0 

(LBS) 

= 

10804 


= 

0 


13221 

00895 

86258 

93778 

00000 

90000 

23196 

42303 

92440 

11294 

98523 

98803 

05302 

00126774 

60434 

00000 

72565 

00000 

27435 

05403 

00000 

00000 

00000 

00000 

39673 

00000 
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TRIM PARAMETERS 


TRIM PITCH AXIS PARAMETER 

= 

-0.79364 

TRIM ROLL AXIS PARAMETER 

= 

0.00000 

TRIM YAW AXIS PARAMETER 

= 

0.00000 

TRIM THRUST PARAMETER 

= 

0.22509 

CONTROL VARIABLES 

ELEVATOR 

= 

0.06377 

THROTTLE 


0.22509 

SPEED BRAKE 

= 

0.00000 


OBSERVATION VARIABLES 



AN 

= 

0.98522771 

AY 

SB 

0.00000000 
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NON-DIMENSIONAL STABILITY AND CONTROL DERIVATIVES FOR CASE #2 
LINEARIZER TEST AND DEMONSTRATION CASES 


w 
w u 
q 2 


lO O O O O O 

*- 0 00 0 0 

I + + + + + 

Q Q Q Q Q Q 

1^-0 0 0 0 0 

CO O O O O O 

P-OOOOO 

m o o o o o 

T- o o o o o 

• ••••• 


o *- o o o 

00000 
+ \ + + + 

q q a Q a 

00000 
o ro o o o 

00000 
0^000 
o r- o o o 

• • • • • 


■^*0000 


000000 


o o 
o o 
+ + 
Q Q 
o o 
o o 
o o 
o o 
o o 


000 


Eh 

2 

H 


*-o*-oor^ooo*-o 
0000*— 000000 
1 + + +1 I+ + + + + 
qqoqqqqqqqq 
oooot-oooooo 

OOCNO-sTLnOOOfNO 

coOroOCNOOOOroO 

t^OCNOO^l^OOCNO 

mor^o<tfi^aooor^o 

«•••••••«•* 

*— Or— ooo^oo*— o 


e> 

s 


a 


CN O 

0 o 

1 + 

Q Q 
O O 
O O 
P- O 
00 O 
O O 


o o 
o o 


ooooocooooo 


CN O CN 
OOO 

I + I 

a a a 

000 

r O ^ 

cn o co 

CO o o 

o o ^ 

• • • 

■<r o ko 


O EHj 
2 
H 


VO 

*- O 
I I 

Q Q 
m o 
co v- 
m cn 
cn r- 

CN CO 

• • 

*- CO 
I 


o o 
o o 
+ + 
a a 
o o 
o o 
o o 
o o 
o o 


o 

o 

+ 

a 

o 

o 

o 

o 

o 


0 ^ 000 t -000 


o 

o 

+ 

Q 

o 

o 

o 

o 

o 


o 

o 

+ 

Q 

o 

o 

o 

o 

o 


000 


w 

.-q 

CQ 

< 

Eh 

CO 

a 

2 

o 

2 

O 


O 

2 

M 

2 


O W 
^ 2 


cNOOOor- 
o o o o *- o 

I + + + I I 

Q Q Q O O Q 

OOOOOOP- 

•^r o co 0 o 

o o m o m m 

cn o o o co cn 

CN O CO O ^ P 

• ••••• 

I 1 o co O 10 10 

I I 


r- O O *- O 

00000 

1 + + + + 

Q Q Q Q Q 

OOOOO 
CN O O P' O 

00 O O 00 O 

CO O O 00 O 

O O O *- O 

• • • • • 

r- O O *- O 

I I 


*- O r- 
OOO 

Q S Q 

OOO 
00 o o 
cn o m 
m o p 
o o *- 
• • • 
O O N* 
I I 


u 

2 

M 

►Q 

2 


vc*— O*— OOOt— OOO 

r -0000000000 

I I + I + + + t + + + 
QQOQQQQQQQQ 
r^oooooooooo 

I^OOOOOOLTIOOO 

moooooo^ooo 

unoooooof’iooo 

(NOOtnOOOrooOO 

• •••••••••• 

*— CNO*— OOO*— OOO 


o 

o 

+ 

Q 

o 

o 

o 

o 

o 


OOO 


o 


2 

>H 

a 

o 

2 

w 

< 

2 

< 

w 

2 


dP 

IT) 

CO 


co 2 
h O 


Eh 

O 

W 

O 

O 

a: 

2 


w 

2 co 
Eh - 


OS W 
O CO 
Cn D 











x-s 




2 

Eh 

u 

O 

CJ 

0 


q 

q 

Eh 

0 

0 




M 

H 

w 

w 

w 

w 


a 

rtj 

Cn 

W 

w 




O 

O 

co 

CO 

CO 

CO 


os 

C8 

— • 

CO 

CO 




2 

2 

CO \ 


\ 

\ 


— - 



\ 

\ 





O 

Eh O 
2 <C 

§ 


Eh 

Cn 





§9 




2 

U 

w 2 

PS 







OS 

BS 




O 

Eh 

W 








— ' 

- — - 




H 

2 

0 













Eh 

O 

M 




OS 








W 

< 

w 

tii 

w 



w 









Eh 


Cii W 

Eh 



00 




Eh 




< 

CO 

2 

W Eh 

kC 

w 

>H 

2 



M 

O 


2 

DQ 

2 



8 g 

2 

Eh 

Eh 

H 

D 

2 



O 

O 

O 

6-* 

O 

O 

Eh 

Eh 

CQ 

W 

J 

CO 

M 


2 

2 

O 


<c 


Eh 

< 

Q 

< 

Eh 

Q 

O 

2 

O 

0 


O 

2 

2 

< 

H 

2 

< 

> 

O 

W 

H 

Eh 

2 J 

Eh 

s 


O 

2 

Eh 

Eh 

2 

Ei 

W 

2 

w 

2 


w 0 

H 

< 

w 

<C 


W 


2 

W 

2 

2 

2 

W 

Eh 

tS3 PS 

2 

>H 

> 

2 

< 

2 

< 

<c 

ca 

W 

Eh 

CO 

> 

< 


104 



MATRIX A USING THE FORMULATION OF THE STATE EQUATION: 


TIME 

DERIVATIVES 
WRT : 

ALPHA 

Q 

THETA 

VEL 


TIME 

DERIVATIVES 

WRT: 

ALPHA 

Q 

THETA 

VEL 


DX/DT = A * X + B * U + D * V 
FOR CASE # 2 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 

ALPHA Q THETA VEL 


-0. 120900D+01 
-0.149189D+01 
0.000000D+00 
-0 . 576868D+02 


0. 1 00000D+0 1 
-0.221451D+01 
0. 100000D+01 
0.000000D+00 


-0.575730D-02 
0. 189640D-01 
0 . 000000D+00 
-0.316251D+02 


-0 . 70 1975D-04 
0.231368D-03 
0 . 000000D+00 
-0.460435D-02 


MATRIX B USING THE FORMULATION OF THE STATE EQUATION: 
DX/DT = A * X + B * U + D * V 
FOR CASE # 2 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 

ELEVATOR THROTTLE SPEED BRAKE 


-0. 141961D+00 
-0 . 220778D+02 
O.OOOOOOD+OQ 
-0. 105186D+02 


0 . 448742D-03 
-0. 147812D-02 
0. 000000 D+ 00 
0 . 343 1 62D+02 


-0.928932D-02 
-0. 135074D+02 
0 . 000000D+00 
-0 . 1 55832D+02 



MATRIX D USING THE FORMULATION OF THE STATE EQUATION: 
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MATRIX H USING THE FORMULATION OF THE OBSERVATION EQUATION: 


Y = H * X + F * U + E * V 


FOR CASE # 2 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 



ALPHA 

Q 

THETA 

VEL 

AN 

0 . 350424D+02 

0.277556D-16 

-0.632324D-02 

0 . 203434D-02 

AY 

0 . 000000D+00 

0 . 000000D+00 

0.000000D+00 

0.2 1 2306D-16 


MATRIX F USING THE 

FORMULATION OF 

THE OBSERVATION 

EQUATION : 


Y = 

H * X + F * 

U + E * V 



AN 

AY 


FOR CASE # 2 

LINEARIZER TEST AND DEMONSTRATION CASES 


FOR THE PROJECT: 
USER'S GUIDE 


ELEVATOR 

0.41 1 323D+0 1 
0 . 000000D+00 


THROTTLE 

0.492845D-03 
0 . 000000D+0Q 


SPEED BRAKE 

0 . 263288D+00 
0 . 000000D+00 
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MATRIX E USING THE FORMULATION OF THE OBSERVATION EQUATION: 
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TRIM CONDITIONS FOR CASE # 1 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 

LEVEL TURN WHILE VARYING ALPHA 
TRIM ACHIEVED 


COEFFICIENT OF LIFT 


= 

0.40144 

COEFFICIENT OF DRAG 


= 

0.03058 

LIFT 

(LBS) 

= 

134741 .68807 

DRAG 

(LBS) 

= 

10265.70661 

ALTITUDE 

(FT) 

= 

20000.00000 

MACH 


= 

0.90000 

VELOCITY 

(FT/SEC) 

= 

933.23196 

EQUIVALENT AIRSPEED 

(KTS) 

= 

403.42303 

SPEED OF SOUND 

(FT/SEC) 

= 

1036.92440 

GRAVITATIONAL ACCEL 

(FT/SEC* *2) 

= 

32.11294 

NORMAL ACCELERATION 

(G-S) 

= 

3.00163 

LOAD FACTOR 


= 

2.99995 

DYNAMIC PRESSURE 

(LBS/FT**2 ) 

= 

552.05302 

DENSITY 

( SLUG/FT**3 ) 

= 

0.00126’ 

WEIGHT ( ©ALTITUDE ) 

(LBS) 

= 

44914.60434 

BETA 

(DEG) 

= 

0.03193 

ALPHA 

(DEG) 

= 

2.66824 

PHI 

(DEG) 

= 

70.62122 

THETA 

(DEG) 

= 

0.91607 

ALTITUDE RATE 

(FT/ SEC) 

= 

0.00000 

GAMMA 

(DEG) 

= 

0.00000 

ROLL RATE 

(DEG/SEC) 

= 

-0.08951 

PITCH RATE 

(DEG/SEC) 

= 

5.28086 

YAW RATE 

(DEG/SEC) 

= 

1.85749 

THRUST 

(LBS) 

= 

10277.03515 

SUM OF THE SQUARES 


= 

0.00000 
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TRIM PARAMETERS 


TRIM PITCH AXIS PARAMETER 
TRIM ROLL AXIS PARAMETER 
TRIM YAW AXIS PARAMETER 
TRIM THRUST PARAMETER 

CONTROL VARIABLES 

ELEVATOR 

THROTTLE 

SPEED BRAKE (DEG) 

OBSERVATION VARIABLES 

AN 

AY 


-0.66958 

-0.01526 

-0.02125 

0.21410 


0.05380 

0.21410 

0.00000 


3.00163339 GS 
.94136286 GS 


VEHICLE STATIC MARGIN IS 3.5% MEAN AERODYNAMIC CHORD STABLE AT THIS 
FLIGHT CONDITION 
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TRIM CONDITIONS FOR CASE # 2 

LINEARIZER TEST AND DEMONSTRATION CASES 

FOR THE PROJECT: 

USER'S GUIDE 

STRAIGHT AND LEVEL TRIM WHILE VARYING ALPHA 


TRIM ACHIEVED 


COEFFICIENT OF LIFT 
COEFFICIENT OF DRAG 
LIFT 
DRAG 

ALTITUDE 

MACH 

VELOCITY 

EQUIVALENT AIRSPEED 
SPEED OF SOUND 
GRAVITATIONAL ACCEL 
NORMAL ACCELERATION 
LOAD FACTOR 
DYNAMIC PRESSURE 
DENSITY 

WEIGHT (©ALTITUDE) 

BETA 

ALPHA 

PHI 

THETA 

ALTITUDE RATE 
GAMMA 
ROLL RATE 
PITCH RATE 
YAW RATE 
THRUST 

SUM OF THE SQUARES 


(LBS) 

(LBS) 

(FT) 

(FT/SEC) 

(KTS ) 

(FT/ SEC) 

(FT/ SEC* *2) = 

(G-S) 

(LB/FT**2) 
(SLUG/FT* *3) = 
(LBS) 

(DEG) 

(DEG) 

(DEG) 

(DEG) 

(FT/ SEC) 

(DEG) 

(DEG/SEC) 

(DEG/SEC) 

(DEG/SEC) 

(LBS) 


0.13221 

0.00895 

44376.86258 

3004.93778 

20000.00000 

0.90000 

933.23196 

403.42303 

1036.92440 

32.11294 

0.98523 

0.98803 

552.05302 

0.00126774 

44914.60434 

0.00000 

-0.72565 

0.00000 

9.27435 

162.05403 

10.00000 

0.00000 

0.00000 

0.00000 

10804.39673 

0.00000 



TRIM PARAMETERS 


TRIM PITCH AXIS PARAMETER 

= 

-0.79364 

TRIM ROLL AXIS PARAMETER 

= 

0.00000 

TRIM YAW AXIS PARAMETER 

= 

0.00000 

TRIM THRUST PARAMETER 

= 

0.22509 

CONTROL VARIABLES 

ELEVATOR 

= 

0.06377 

THROTTLE 

= 

0.22509 

SPEED BRAKE 

= 

0.00000 


OBSERVATION VARIABLES 



AN 

= 

0.98522771 

AY 

= 

0.00000000 


VEHICLE STATIC MARGIN IS 3.5% MEAN AERODYNAMIC CHORD STABLE AT 
THIS FLIGHT CONDITION. 




APPENDIX J: EXAMPLE USER-SUPPLIED SUBROUTINES 


The following subroutines are examples of the user-supplied routines that provide 
the aerodynamic, control, engine, and mass and geometry models to LINEAR. These sub- 
routines are based on an F-15 simulation and are typical of the routines needed to 
interface LINEAR to a set of nonlinear simulation models. These subroutines are 
meant to show the use of the named common blocks to communicate between LINEAR and 
the user's routines. These subroutines are used with all examples in this report. 
Included with this report are microfiche listings of these subroutines. 


Aerodynamic Model Subroutines 

The following two subroutines define a linear aerodynamic model. Even though 
this model is greatly simplified from the typical nonlinear aerodynamic model, the 
example shows the functions of the subroutines ADATIN and CCALC. 


SUBROUTINE ADATIN 


C 

C EXAMPLE AERODYNAMIC DATA DEFINITION (OR INPUT) SUBROUTINE 

C 

C.... ROUTINE TO DEFINE STABILITY AND CONTROL DERIVATIVES FOR THE 
C AERODYNAMIC MODEL 

C 
C 

C COMMON BLOCK TO COMMUNICATE AERODYNAMIC DATA BETWEEN THE 

C SUBROUTINES ADATIN AND CCALC 

C 
C 


COMMON /ARODAT/ CLB 

,CLP 

,CLR 

/ 

CLDA 

,CLDR 

,CLDT 

/ 

CMO 

,CMA 

/CM Q 

,CMAD , 

CMDE 

,CMSB 

/ 


CNB 

,CNP 

,CNR 

# 

CNDA 

,CNDR 

,CNDT 

/ 

CDO 

,CDA 

, CDDE 

, CDSB , 

CLFTO 

,CLFTA 

,CLFTQ 

,CLFTAD, 

CLFTDE 

,CLFTSB 

/ 


CYB 

,CYDA 

,CYDR 

,CYDT 


C 

C . . . . ROLLING MOMENT DERIVATIVES . 

C 

C STABILITY DERIVATIVES WITH RESPECT TO 

C SIDESLIP, ROLL RATE, YAW RATE . . . 

C 



CLB 

= -1.3345 

D-0 1 


CLP 

= -2.0000 

D-0 1 

C 

CLR 

= 1.5099 

D-0 1 

C 


CONTROL 

DERIVATIVES WITH RESPECT TO 

C 

C 


AILERON , 

RUDDER, DIFFERENTIAL TAIL 
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noon oooooo ooon 


CLDA = 2.6356 D-02 

CLDR = -2.3859 D-03 
CLDT = 4.0107 D-02 

C 

C.... PITCHING MOMENT DERIVATIVES. 

C 

C STABILITY DERIVATIVES WITH RESPECT TO ANGLE OF 

C ATTACK, PITCH RATE, ANGLE OF ATTACK RATE . . . 

C 

CMO = 4.2204 D-02 

C 

CMA = -1.6882 D-0 1 
CMQ = 3.8953 

CMAD = -1 . 1887 D+0 1 

CONTROL DERIVATIVES WITH RESPECT TO 
ELEVATOR, SPEED BRAKE . . . 

CMDE = -6.9528 D-01 
CMSB = -4.1750 D-01 

YAWING MOMENT DERIVATIVES. 

STABILITY DERIVATIVES WITH RESPECT TO 
SIDESLIP, ROLL RATE, YAW RATE . . . 


CNB 

= 1.2996 

D-01 

CNP 

= -3.3721 

D-02 

CNR 

= -4.0471 

D-01 


CONTROL 

DERIVATIVES WITH RESPECT TO 


AILERON, 

RUDDER, DIFFERENTIAL TAIL 

CNDA 

= 2.1917 

D-03 

CNDR 

= -6.9763 

D-02 

CNDT 

= 3.0531 

D-02 


C 

C.... COEFFICIENT OF DRAG DERIVATIVES. 

C 

C STABILITY DERIVATIVES WITH RESPECT TO 

C ANGLE OF ATTACK . . . 

C 

CDO = 1.0876 D-02 

C 

CDA = 3.7257 D-01 

C 
C 

C CONTROL DERIVATIVES WITH RESPECT TO 

C ELEVATOR, SPEED BRAKE . . . 

C 

CDDE = 4.3831 D-02 

CDSB = 6.4935 D-02 

C 
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C.... COEFFICIENT OF LIFT DERIVATIVES. 

C 

C STABILITY DERIVATIVES WITH RESPECT TO ANGLE OF 

C ATTACK, PITCH RATE, ANGLE OF ATTACK RATE . . . 

C 

CLFTO = 1.5736 D-01 

C 

CLFTA = 4.8706 

CLFTQ = -1.7232 D+01 
CLFTAD= 1.7232 D+01 
C 

C CONTROL DERIVATIVES WITH RESPECT TO 

ELEVATOR, SPEED BRAKE . . . 

CLFTDE= 5.7296 D-01 
CLFTSB= 3.7492 D-02 


.... SIDEFORCE COEFFICIENT DERIVATIVES. 

STABILITY DERIVATIVES WITH RESPECT TO 
SIDESLIP . . . 

CYB = -9.7403 D-01 

CONTROL DERIVATIVES WITH RESPECT TO 
AILERON, RUDDER, DIFFENTIAL TAIL 
C 

CYDA = -1.1516 D-03 
CYDR = -1.5041 D-01 
CYDT = -7.9315 D-02 
C 

RETURN 

END 

SUBROUTINE CCALC 
C 

C.... EXAMPLE AERODYNAMIC MODEL. 

C 

C . . . . ROUTINE TO CALCULATE THE AERODYNAMIC FORCE AND MOMENT COEFFICIENTS . 
C 
C 
C 
C 
C 

COMMON 
COMMON 
C 

COMMON 
COMMON 

COMMON 
COMMON 
C 


COMMON BLOCKS CONTAINING STATE, CONTROL, AND AIR 
DATA PARAMETERS 


/DRVOUT/ F( 13 ) ,DF( 13) 

/CONTRL/ DC (30) 

/DATAIN/ S, B ,CBAR, AMSS , AIX, AIY , AIZ , AIXZ , AIXY , AIYZ , AIXE 
/TRIGFN/ SINALP , COSALP , SINBTA , COSBTA , SINPHI , COSPHI , 
SINPSI , COSPSI , SINTHA,COSTHA 
/SIMOUT/ AMCH , QBAR, GMA , DEL , UB , VB , WB , VEAS , VCAS 
/CGSHFT/ DELX , DELY , DELZ 
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C COMMON BLOCK TO OUTPUT AERODYNAMIC FORCE AND MOMENT 

COEFFICIENTS 


COMMON / CLCOUT/ CL ,CM ,CN ,CD ,CLFT ,CY 


COMMON BLOCK TO COMMUNICATE AERODYNAMIC DATA BETWEEN 
THE SUBROUTINES ADATIN AND CCALC 


COMMON /ARODAT/ CLB 

,CLP 

f CLR 

/ 

CLDA 

,CLDR 

, CLDT 

/ 

CMO 

,CMA 

,CMQ 

,CMAD , 

CMDE 

, CMSB 



CNB 

,CNP 

,CNR 

t 

CNDA 

,CNDR 

,CNDT 

/ 

CDO 

,CDA 

,CDDE 

z CDSB , 

CLFTO 

,CLFTA 

,CLFTQ 

/CLFTAD, 

CLFTDE 

/CLFTSB 



CYB 

,CYDA 

,CYDR 

,CYDT 


EQUIVALENCE VARIABLE NAMES 


EQUIVALENCE 

(T 

, F ( 

D), 




• 

(P 

, F( 

2) ),(Q 

, F( 3)), 

(R 

, F ( 4)) 


(V 

/ F( 

5)), (ALP 

, F( 6)), 

( BTA 

, F( 7)) 


(THA 

, F ( 

8) ) , (PSI 

, F( 9)), 

(PHI 

, F ( 10) ) 


(H 

, F( 1 1 ) ) » (X 

, F( 1 2 ) ) , 

(Y 

, F ( 13) ) 


(TDOT 

,DF( 

D), 





( PDOT 

,DF( 

2 ) ) , ( QDOT 

,DF( 3)), 

(ROOT 

,DF( 4)) 


(VDOT 

,DF( 

5 ) ) , ( ALPDOT 

,DF( 6)), 

( BTADOT 

,DF( 7)) 


( THADOT 

,DF( 

8) ), (PSIDOT 

,DF( 9)), 

(PHIDOT 

,DF( 10) ) 


(HDOT 

,DF( 11)) , (XDOT 

,DF( 12) ), 

( YDOT 

,DF( 13) ) 

EQUIVALENCE 

(DA 

,DC( 

1 ) ) / ( DE 

,DC( 5)), 

(DT 

DC( 8)), 

• 

(DR 

,DC( 

9) ) , (DSB 

,DC( 10) ) 



COMPUTE TERMS NEEDED 

WITH 

ROTATIONAL 

DERIVATIVES. 



V2 =2.0 *V 

B2V =B /V2 

C2V =CBAR /V2 

. ... ROLLING MOMENT COEFFICIENT. 


CL 

=CLB 

* BTA 

+CLDA 

*DA 

+CLDR 

*DR 

+CLDT 


+B2V 

* ( CLP 

*P 

+CLR 

*R 

) 


PITCHING MOMENT COEFFICIENT. 





CM 

=CMO 

+ CMA 

*ALP 

+CMDE 

*DE 

+CMSB 

*DSB 


+C2V 

* (CMQ 

*Q 

+CMAD 

* ALPDOT) 




C 

C.... YAWING MOMENT COEFFICIENT. 
C 


*DT 
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CN =CNB * BTA +CNDA *DA +CNDR *DR +CNDT *DT 

. +B2V * ( CNP *P +CNR *R ) 

C 

C.... COEFFICIENT OF DRAG. 

C 

CD =CDO + CDA *ALP +CDDE *DE +CDSB *DSB 

C 

C.... COEFFICIENT OF LIFT. 

C 

CLFT =CLFTO + CLFTA *ALP +CLFTDE*DE +CLFTSB*DSB 

. +C2V * ( CLFTQ *Q +CLFTAD * ALPDOT ) 

C 

C SIDEFORCE COEFFICIENT. 

C 

CY =CYB *BTA +CYDA *DA +CYDR *DR +CYDT *DT 

C 

RETURN 

END 


Control Model Subroutine 

The following subroutine UCNTRL models the gearing between the control inputs 
(that is, trim parameters 6a, 6e, 6r, and thr) and the control surfaces. This pro- 
vides LINEAR with a generic interface capable of handling most airborne vehicles. 
Figure 8 shows the control system implemented in this example subroutine UCNTRL. 


DES 

DAS 

DRS 

THRSTX 



UCNTRL 


Figure 8 . Gearing model in example 
UCNTRL subroutine . 
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SUBROUTINE UCNTRL 


C 

C 

C 

C 

C 

C 

C 

C 


EXAMPLE TRIM/CONTROL SURFACE INTERFACE ROUTINE. 

ROUTINE TO CONVERT TRIM INPUTS INTO CONTROL SURFACE DEFLECTIONS. 

INPUT COMMON BLOCK CONTAINING TRIM PARAMETERS 
COMMON / CTPARM/ DES ,DAS ,DRS ,THRSTX 


C 

C 

C 

C 

C 

C 


C 

C 

c . . . . 

c 


c 


OUTPUT COMMON BLOCK CONTAINING CONTROL SURFACE DEFLECTIONS. 
COMMON /CONTRL/ DC (30) 

EQUIVALENCE VARIABLE NAMES. 

EQUIVALENCE (DA ,DC( 1 ) ) , ( DE ,DC( 5)),(DT ,DC(8)), 

(DR ,DC( 9 ) ) , ( DSB ,DC( 10 ) ) , (THR ,DC(12)) 

DATA DGR /57. 29578/ 

CONVERT FROM INCHES OF STICK AND PEDAL TO DEGREES OF SURFACE. 


DA 

=DAS 

*( 20.0 

/ 

4.0 ) 

DE 

=DES 

*(-25.0 

/ 

5.43) 

DR 

=DRS 

*( 30.0 

/ 

3.25) 


SET DIFFERENTIAL TAIL BASED ON AILERON COMMAND. 
DT =DA / 4.0 


.... CONVERT THRUST TRIM PARAMETER TO PERCENT THROTTLE COMMAND. 
THR =0.0 

IF(THRSTX.GE. 0.0 ) THR =THRSTX 
C 

C.... USE SPEED BRAKE IF NEEDED. 

C 

DSB =0.0 

IF (THRSTX.LT. 0.0 ) DSB =THRSTX* ( -45 . 0 ) 

C 

C . . . . CONVERT SURFACE COMMANDS TO RADIANS . 


DA 

=DA 

/DGR 

DE 

=DE 

/DGR 

DR 

=DR 

/DGR 

DT 

=DT 

/DGR 

DSB 

=DSB 

/DGR 


C 

RETURN 

END 
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Engine Model Interface Subroutine 


The following subroutine IFENGN both defines an engine model and provides the 
interface to LINEAR. In normal usage, this subroutine would call subroutines that do 
the detailed thrust, engine rotation, and fuel consumption modeling; the information 
from these subroutines would be transferred into the /ENGSTF/ common block. 


C 

C . . . . 

c 

c . . . . 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c . . . . 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE IFENGN 

EXAMPLE SUBROUTINE TO PROVIDE PROPULSION SYSTEM MODEL. 

ROUTINE TO COMPUTE PROPULSION SYSTEM INFORMATION FOR LINEAR 

THIS SUBROUTINE IS THE INTERFACE BETWEEN THE DETAILED ENGINE 
MODELING SUBROUTINES AND LINEAR. 

INPUT COMMON BLOCK CONTAINING INFORMATION ON THRUST 
REQUEST TO ENGINE 

COMMON /CONTRL/ DC(30) 

OUTPUT COMMON BLOCK CONTAINING DETAILED INFORMATION 
ON EACH OF UP TO FOUR SEPARATE ENGINES 

ROUTINE TO COMPUTE ENGINE PARAMETERS: 

THRUST ( I ) = THRUST CREATED BY EACH ENGINE 
TLOCAT( I , J) = LOCATION OF EACH ENGINE IN THE X-Y-Z PLANE 
XYANGL(I) = ANGLE IN X-Y BODY AXIS PLANE AT WHICH EACH 
ENGINE IS MOUNTED 

XZANGL(I) = ANGLE IN X-Z BODY AXIS PLANE AT WHICH EACH 
ENGINE IS MOUNTED 

TVANXY ( I ) = ANGLE IN THE X-Y ENGINE AXIS PLANE OF THE 
THRUST VECTOR 

TVANXZ(I) = ANGLE IN THE X-Z ENGINE AXIS PLANE OF THE 
THRUST VECTOR 

DXTHRS ( I ) = DISTANCE BETWEEN THE ENGINE C.G. AND THE 
THRUST POINT 

EIX(I) = ROTATIONAL INERTIA OF EACH ENGINE 
AMSENG ( I ) = MASS OF EACH ENGINE 
ENGOMG ( I ) = ROTATIONAL VELOCITY (RAD/SEC) 

COMMON /ENGSTF/ THRUST ( 4 ) , TLOCAT (4,3), XYANGL ( 4 ) , XZ ANGL ( 4 ) , 

. TVANXY ( 4 ) , TVANXZ ( 4 ) , DXTHRS ( 4 ) , 

. EIX(4 ) ,AMSENG( 4),ENGOMG(4) 



C EQUIVALENCE VARIABLE NAMES 

C 

EQUIVALENCE (THR ,DC(12)) 

C 

C ASSUME THRUST PER ENGINE IS HALF VEHICLE WEIGHT. 

C 

THRUST ( 1 )= 24000.0 *THR 
THRUST (2 )= 24000.0 *THR 
C 

C . . . . LET ALL OTHER PARAMETERS DEFAULT TO ZERO 
C 

RETURN 

END 


Mass and Geometry Model 

The subroutine MASGEO enables the user to model mass, inertia, and center-of- 
gravity changes due to vehicle configuration or mass changes. This particular ver- 
sion is simply a RETURN and END statement, because the aircraft that it models is a 
fixed-wing vehicle, and any mass changes can be done through the input files or even 
interactively. However, this subroutine does provide an elegant way of including a 
model of those changes for other aircraft. 


C 

C. 

C 

C 


C 


SUBROUTINE MASGEO 


. SUBROUTINE TO COMPUTE THE MASS AND GEOMETRY PROPERTIES OF THE 
AIRCRAFT. 


COMMON /CONPOS/ DAP, FLAT, DATRIM,DEP, FLON , DETRIM, DRP, FPRF ,DRTRIM, 
. DSBP , DFP , PLAPL , PLAPR, THSL , THSR 

COMMON /CONTRL/ DC (30) 

COMMON /DATAIN/ S , B , CBAR, AMSS , AIX , AIY , AIZ , AIXZ , AIXY , AIYZ , AIXE 
COMMON /SIMOUT/ AMCH ,QBAR, GMA , DEL, UB , VB , WB , VEAS , VCAS 
COMMON /DRVOUT/ F ( 1 3 ) , DF ( 1 3 ) 


EQUIVALENCE 


(T 

(P 

(V 

(THA 

(H 

(TDOT 
( PDOT 
(VDOT 


F ( 
F( 
F( 
F( 


1 )># 

2 ) ) , (Q 
5 ) ) , ( ALP 
8) ) , (PSI 


, F( 1 1 ) ) , (X 
, DF( 1)), 

,DF( 2 ) ) , ( QDOT 
,DF( 


F( 3 ) ) , ( R 
F( 6 ) ) , ( BTA 
F( 9)), (PHI 
F( 1 2 ) ) , ( Y 


F( 4)), 
F( 7)), 
F ( 10) ) , 
F(13)), 


, DF( 3 ) ) , ( RDOT ,DF( 4)), 
5 ) ) , ( ALPDOT , DF ( 6 ) ) , ( BTADOT , DF ( 7)), 
( THADOT , DF ( 8 ) ) , ( PSIDOT ,DF( 9 ) ) , ( PHIDOT,DF( 1 0 ) ) , 
(HDOT , DF( 11)), (XDOT , DF ( 1 2 ) ) , ( YDOT ,DF(13)) 


RETURN 

END 
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